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Large-distance and long-time asymptotic behavior of the reduced 
density matrix in the non-linear Schrodinger model. 



o 

(N 



in 



K. K. Kozlowskfl 



Abstract 



Starting from the form factor expansion in finite volume, we derive the multidimen- 
sional generalization of the so-called Natte series for the time and distance dependent 
reduced density matrix at zero-temperature in the non-linear Schrodinger model. This 
representation allows one to read-off straightforwardly the long-time/large-distance 
asymptotic behavior of this correlator. This method of analysis reduces the complexity 
of the computation of the asymptotic behavior of correlation functions in the so-called 
interacting integrable models, to the one appearing in free fermion equivalent models. 
We compute explicitly the first few terms appearing in the asymptotic expansion. Part 
of these terms stems from excitations lying away from the Fermi boundary, and hence 
, go beyond what can be obtained by using the CFT/Luttinger liquid based predictions. 
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5 1 Introduction 

^ One-dimensional quantum models with a gapless spectrum are believed to be critical at zero temperature. In 
• ^ other words, in these models, the ground state expectation values of products of local operators should decay, 
^ for large distances of separation between the operators, as some power-law in the distance. It is also believed 
d ' that, for a generic class of Hamiltonians, the actual value of the exponents governing this power-law decay, the 
so-called critical exponents, does not depend on the microscopic details of the interactions in the model, but only 
on its overall symmetries [21^ 2A1- Therefore two models belonging to the same universality class should be 
characterized by the same critical exponents. 

It has been argued that the equal-time correlation functions in quantum critical one-dimensional models exhibit 
conformal invariance in the large-distance regime |74|. Hence, its appears plausible to infer their large-distance 
asymptotics from those of the associated conformal field theory (CFT). The central charge of the CFT lying in the 
universality class of the model can be deduced from the finite-size corrections to the ground state energy lISlfTTI. 
The possibility to compute such finite-size corrections for many integrable models allowed the identification of the 
central charge and scaling dimensions leading to the predictions for the critical exponents E [171 ttH EH EJ l84l 
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of the long-distance asymptotics. We remind that it is also sometimes possible to give predictions for the critical 
exponents |[25l l26l l67l by putting the model in correspondence with a Luttinger liquid ETl . 

Due to their wide applicability and relative simplicity, it is more than desirable to test these CFT/Luttinger 
liquid based predictions versus some exact calculations carried out on such models; this starting from first principle 
and in such a way that no approximation (apart from assuming a large distance) is made to the very end. Such 
exact computations have been carried out in the 70's and 80's on various two-point functions appearing in free 
fermion equivalent models such as the Ising ifTBHTOi rTTl. the XY model at the critical magnetic field Il68ll69l or the 
impenetrable Bose gas |83|. The latter approaches were then made much more systematic (and also simplified) 
with the occurrence of a Riemann-Hilbert based approach to free-fermion models' asymptotics |[30l together with 
the development of the non-linear steepest-descent method |[T6l . Indeed, the latter constitutes a relatively simple 
and systematic tool for carrying out the asymptotic analysis |[29l [311 [32l of Riemann-Hilbert problems associated 
with Fredholm determinants representing the correlators in free fermionic models. 

However, obtaining long-distance asymptotic expansions of two-point functions for models not equivalent to 
free fermions faced several additional problems of technical nature. This fact takes its roots in that even obtaining 
explicit expressions for the coiTclation functions in the so-called interacting integrable models demands to over- 
come new types of combinatorial intricacy that disappears when dealing with free fermion equivalent models. The 
first approach to the problem of computing correlators out of the free fermion point can be attributed to Izergin 
and Korepin |[35l[36ll . These authors managed to construct certain series representations for the correlation func- 
tions of the non-linear Schrodinger model and the XXZ spin- 1/2 chain. However, the n"^ summand appearing 
in these series was only defined implicitly by induction. Low-n calculation allowed them for an effective per- 
turbative characterization of a vicinity of the free fermion point. First manageable expressions for correlators at 
zero temperature in an interacting integrable model were obtained by Jimbo, Miki, Miwa, Nakayashiki through 
the vertex operator approach. They have provided multiple integral representations for the matrix elements of 
the so-called elementary block^ in the massive Ii39i1 regime of the infinite XXZ chain. These results where later 
extended to the massless regime of this chain HOl or to a half-infinite chain subject to a longitudinal magnetic 
field acting on one of its ends Il38l . The multiple integral representations were then reproduced, in the framework 
of the algebraic Bethe Ansatz by Kitanine, Maillet and Terras for the massive and massless regime of the periodic 
XXZ chain [52]. These two series of works opened a way towards a systematic and effective computation of 
various types of multiple integral and/or combinatorial representations for the correlation functions in numerous 
integrable models. In particular, it was possible to derive effective representations in the case of finite temperature 
|[23]| . non-equal times [49], models in finite volume |[37i . higher spin chains [15],... These results should be seen 
as of uttermost importance from the conceptual point of view: the multiple integral representations for the correla- 
tors of interacting integrable models naturally provides an interpretation for these objects as a new class of special 
functions (of the distance, time, coupling constants, ...). However, the complexity of the integrands appearing in 
such multiple integral representations makes the thorough description (computation at certain specific values of 
the distance/coupling or extracting their large-distance/long-time behavior etc..) of these new special functions a 
quite challenging problem. Many investigations that followed where oriented towards a better understanding of 
these special functions. In particular, it was observed that the multiple integral representations for the elementary 
blocks of the XXZ chain can be reduced to one dimensional integrals by a case-by-case analysis |[T0l[47l[50ir75l . 
This observation led to the proof that it is possible to separate the multiple integrals representing the elementary 
blocks of the XXZ chain on the algebraic level [7J. In its turn, this led to the discovery of a Grassmann structure 
in the XXZ chain [(6l[8l[3. Among many other developments such as the possibility to compute the one-point 
functions of the sine-Gordon model ll4n[42i . the existence of such a Grassmann structure constitutes a promising 
direction towards bringing the complexity of the analysis of the correlation function in the XXZ chain to the one 
of a free fermion problem. A completely other method for reducing the complexity of the multiple integral repre- 

' these constitute a basis on wiiicli it is possible to decompose all correlation functions of the model 
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sentations for the correlation functions was the so-called dual field approach |[59l . It led to representations for the 
correlators in terms of expectation values in an infinite dimensional Hilbert space of unbounded operator valued 
Fredholm determinants 1561 l57l . However, apart from convergence issues posed by such an infinite dimensional 
framework, the main problem of that method was posed by the non-commutativity of the dual field's vacuum 
expectation values and the asymptotic expansion of a dual-field valued Fredholm determinant. Its and Slavnov 
|[33l carried out, on a formal level, such a dual field-based analysis for the large-distance/long-time decay of the 
so-called one-particle reduced density matrix at finite temperature in the non-linear Schrodinger model (NLSM). 
They have been able to provide operator valued expressions for the correlation length. The dual field vacuum ex- 
pectation values where computed in |80|, which led to a description of the correlation length in terms of a solution 
to a non-linear integral equation. We would like to mention that until recently, although formal, the dual field 
approach was the only approach alternative to a CFT/Luttinger liquid based correspondence that allowed one to 
write down certain predictions for the critical exponents. 

There have also been developments aiming at obtaining alternative types of effective series of multiple integral 
representations for the distance dependent two-point functions BSl |5TI . The guideline being a construction of 
such a representation that would allow one to carry out a long-distance asymptotic analysis of the two-point 
functions. This project has met a success in [45 1. This article developed a first fairly rigorous method allowing 
one to compute, starting from "first principles", the long-distance asymptotic behavior of the spin-spin correlation 
function in the massless regime of the XXZ spin- 1/2 chain. This method rehed on a few conjectures relative to the 
permutation of symbols, summability of the remainders, convergence of the obtained series representations, but 
was rigorous otherwise. These last results not only confirmed the CFT/Luttinger liquid-based predictions for the 
critical exponents in this model but also provided explicit expressions for the amplitudes in front of the power-law 
which, in their turn, cannot be predicted by universality arguments. These explicit formulae for the amplitudes 
were then identified with certain, properly normalized in respect to the size of the system, form factors of the 
spin operators Il44ll46l . This identification allowed one to point out the universality in the power of the system- 
size that one uses for normalizing the form factor associated with the amplitudes. The aforementioned method 
of asymptotic analysis was applied recently to study the long-distance asymptotic behavior of the correlation 
functions at finite temperature in the NLSM [63]. 

The large-distance/long-time asymptotic behavior of the correlation functions in massless one dimensional 
quantum models goes beyond the predictions stemming from a correspondence with a CFT/Luttinger liquid. 
Hence, this constitutes a clear motivation for obtaining such an asymptotic behavior from exact considerations 
on some integrable model as this could help to understand their structure and origin in the general case where 
exact computations are not feasible. We would like to mention that there already exists several exact results rel- 
ative to this regime of the asymptotics in the case of free fermion equivalent models Il32l l68l |69l 1721 . We also 
would like to remind that there has been proposed recently ||28l a non-Hnear Luttinger liquid theory allowing one 
to predict the leading power-law behavior near the edges of the purely hole or particle specta for dynamic struc- 
ture factors and spectral function^ at low energy and momentum. This approach has been combined with Bethe 
Ansatz considerations to propose values for the amplitudes in front of this behavior [ 12 |. 

This article develops a method allowing one to compute the zero-temperature asymptotic behavior of the 
correlation functions in integrable models starting from the form factor expansion for two-point functions. The 
fact that we build our method on the form factor expansion allows us to include the time-dependence and hence 
access to the large-distance and long-time asymptotic behavior. The method has been introduced recently on 
the example of the current-current correlators [65 1. Here, we provide many elements of rigor to the method and 
treat the example of the one particle reduced density matrix in the non-linear Schrodinger model. We would 
like to stress that this method of asymptotic analysis not only allows us to carry out the analysis in the large- 

^These quantities refer to space and time Fourier transforms of particular two-point functions 
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distance/long-time regime but also constitutes an important technical and computational simplification of the 
approach proposed in Il45l . It has also the advantage of being applicable to a much wider class of algebraic Bethe 
Ansatz solvable models as it solely relies on the universal structure of the form factors in these models. All the 
more than the number of models where these have been determined is constantly growing lfT4ll36ll53ll6Tir731 . The 
main result of this paper can be summarized as follows. We provide a method for constructing a new type of series 
representation for the correlation functions of integrable models, that we call multidimensional Natte series. This 
representation is THE one that is fit for an asymptotic analysis, as the first few terms of the asymptotic expansion 
can be simply read-off without any effort by looking at the terms of the series. Moreover, the computation of the 
higher order asymptotics effectively boils down to the case of a free fermionic model (ie computation of subleading 
asymptotics of the Fredholm determinant of an integrable integral operator) and thus bears the same combinatorial 
complexity. The main implication of our result for physics is that the asymptotics in the time-dependent case are 
not only driven by excitations on the Fermi boundary (the latter coincides with the region of the spectrum that 
can be taken into account by using CFT/Luttinger liquid-based predictions), but also by excitations around the 
saddle-point Aq of the "plane-wave" combination xp (A) - ts {A) of the dressed momentum p and dressed energy 
E of the excitations. Also, we provide explicit expressions and identify the associated amplitudes with the infinite 
volume limit of the properly normalized in the size of the system form factors of the field. 

We stress that although we have been able to set our method in a more rigorous framework then it was done 
in P31I631 . we still have to rely on a few conjectures. More precisely, we have been able to split the asymptotic 
analysis part from the one of proving the convergence of certain series of multiple integrals representing the 
correlators. The part related to asymptotic analysis has been set into a rigorous framework. However, in order 
to raise the results of this asymptotic analysis to the level of the two-point function of interest, we still need to 
assume the convergence of the series of multiple integrals we obtain. 

The main novelty of this method is that it provides a systematic way for carrying out the asymptotic analysis of 
multiple integrals or series thereof whose integrands contain some large-parameter dependent driving term being 
dressed up by coupled functions of the integration variables. We provide a setting that allows one to interpret the 
"coupled" case as some deformation of the "uncoupled" one. This deformation is such that, provided one is able to 
carry out the analysis in the "uncoupled" case (but with a sufficiently rich range of functions involved), one is able 
to deform the "uncoupled" asymptotics back to the "coupled" case of interest. It is in this respect that the analysis 
carried out in this article strongly relies on the results obtained in [62] (where the relevant "uncoupled" series of 
multiple integrals of interest has been analyzed) as well as on the fact that correlation functions of generalized free 
fermionic models (which correspond to the "uncoupled" case) are naturally representable in terms of Fredholm 
determinants [60|. 

This paper is organized as follows. In section |2l we remind the definition and main properties of the model. 
We also introduce all the necessary notations allowing us to present the asymptotic behavior of the reduced density 
matrix. In section [3l we present our result and discuss the strategy of our method. Then, in section |4l we outline 
the main properties of the form factors in the model and write down the form factor series for the reduced density 
matrix. We explain how this series can be re-summed into the so-called multidimensional deformation of the 
Natte series. Once that such a representation is built, thanks to the very properties of the Natte series, it is possible 
to literally read-off the first few terms of the asymptotic expansion. We gather all the auxiliary and technical 
results in several appendices. We discuss the large size-behavior of the form factors of the fields in appendix 
lAl In appendix |BJ we derive finite-size Fredholm minor representations for the form factor based expansions of 
certain two-point functions in generalized free fermion models. In appendix O we prove the existence of the 
thermodynamic {ie infinite volume) limit for certain quantities of interest. We also provide various alternative 
expressions for this limit. In appendix |Dl we develop the theory of functional translation in spaces of holomorphic 
functions. The results established in this appendix constitute the main tools of our analysis. They allow for an 
effective separation of variables in the intermediate steps so that one is able to carry out various re-summations of 
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the formulae by building on the results stemming from the generalized free fermion model studied in appendix IbI 

2 The non-linear Schrodinger model 

2.1 The eigenstates and Bethe equations 

The non-linear Schrodinger model corresponds to the Hamiltonian 

L 

Hnls = j [dy^Hy) dy^ (y) + c ^Hy) ^^(y) O (y) ^(y)-h ^^(y) O (y)] dy . (2.1) 


The model is defined on a circle of length L, so that the canonical Bose fields <1)^ are subject to L-periodic 
boundary conditions. In the following, we will focus on the repulsive regime c > in the presence of a positive 
chemical potential h > 0. The Hamiltonian H^ls commutes with the number of particles operator, and thus can 
be diagonalized independently in every sector with a fixed number of particles N. In each of these sectors, the 
model is equivalent to a N-body gas of bosons subject to d-like. repulsive interactions. The corresponding model 
of interacting bosons was first proposed and studied by Girardeau [22] in the c = +oo case and then introduced 
and solved, through the coordinate Bethe Ansatz, by Lieb and Liniger f66\ in the case of arbitrary c. It is also 
possible to build the eigenstates of the Hamiltonian by means of the algebraic Bethe Ansatz. This was first done 
by Sklyanin fTTl directly in the infinite volume. In the case of finite volume L, as observed by Izergin and Korepin 
|[34l . it is possible to put the continuous model on a lattice in such a way that the standard construction f20l of the 
algebraic Bethe Ansatz holds. At the end of the computations, it is then possible to send the lattice spacing to zero 
and recover the spectrum and eigenstates of the continuous model. The fact that this manipulation is indeed fully 
rigorous has been shown by Dorlas |[T9l . 

In the algebraic Bethe Ansatz approach, the Hamiltonian (12.11 ) appears as a member of a one-parameter com- 
muting family of operators A t-^ T (A). It is sometimes useful to consider a yS-deformation of this family TJj (A), 
such that T^(/l)|y5=o = T' i'^)- The common eigenstates | lA/j ({//})) of TpiA) in the A/^^^ -particle sector are param- 
eterized by a set of real numbers which are the unique solutions to the yS-deformed logarithmic Bethe 
equations B1I851 

Lpo (^4) + 2 ^ (/^^" - /^^*) " - J + with po (A) = A and 6 (A) - / In ( ■ (2-2) 

po is called the bare momentum and 9 the bare phase. The set of solutions corresponding to all choices of integers 
e Z such that £1 < • • • < €1^^ yield the complete set of eigenstates in the A^^ -particle sector [,19. 1 . 

In each sector with a fixed number of particles N^, the so-called ground state's Bethe roots are given by the 
solution to (12.21 ) corresponding to the choice of A'^^ consecutive integers = a, with a = I, . . .,Nk and p -0. The 
number corresponding to the number of particles in the overall ground state of Ht^i^ is imposed by the chemical 
potential h and scales with L. It will be denoted by N in the following. One shows that in the thermodynamic limit 
(A'^, L +00 so that N/L — > D) the parameters {Aj}^ associated to this ground state condensate on a symmetric 
interval [-g',t]] called the Fermi zone. 

All other choices of sets of integers {a lead to (yS-deformed) excited states. In principle, these excited states 
can also be found in sectors with a different number N of particles. It is convenient to describe the excited 
states in the language of particle-hole excitations above the A'^*^ -particle /3-deformed ground stat^. Namely, such 

' the yS-deformed ground state corresponds to the choice = a, with a = I, . . . ,Nk 
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an excited state corresponds to a choice of integers €j in (12.21) such that 

€j = j iov j ell; Ni,J\hi,...,h„ and ^h, = Pa for a = 1, . . . , « . (2.3) 

The integers pa and ha are such that pa i \ l\ T^k\ = U, ■ ■ ■ , N,^} and ha e E 1 ; A'^k 3- There is thus a one-to-one 
con^espondence between integers ij and the integers ha and pa describing particle-hole excitations. 

In this picture, the integers ha correspond to holes in the increasing sequence of integers defining the yS- 
deformed ground state roots, whereas pa correspond to extra integers appearing in the equation and can be seen 
as defining some new position of "particles". Given a solution {ncj^" corresponding to a fixed choice of integers 
£i < ■ ■ ■ < it is convenient to introduce their counting function: 

By construction, it is such that ^[[^] (pip) - fa/L, for a = 1, . . . , A'^^. Actually, ^[^i (oj) define^a set of background 
parameters {fia}, a e Z, as the uniquqj solutions to ^[(^ {^a) = ajL. The latter allows one to define the rapidities 
//p,, resp. ^h„, of the particles, resp. holes, entering in the description of jju^^)^'. 

2.2 The thermodynamic Umit 

When the thermodynamic limit of the model is considered, it is possible to provide a slightly more precise de- 
scription of the solution to the Bethe equations for the ground state {/^al^^j as well as for any particle-hole type 
yS-deformed excited states {/^4)^^j above it with N^-N being fixed and not depending on L or N. Introducing the 
counting function for the ground state 

^ ?(w I [Aaf,] = + ^ Z ^ - + ' ^^^"^ = I ' (2-^^ 

a=l 

it can be shown that, in the thermodynamic limit, it behaves as 

Jica) ^^{oj) + (L"M where ^(w) = + ^ and N/L D . (2.6) 

^ ' 2n 2 

There, the O is uniform and holomorphic in a> belonging to a strip of some fixed width around the real axis, 
p is the so-called dressed momentum, defined as the unique solution to the integral equation 

p(A)- j e(A-iu)p'(ju)^= po (A) . (2.7) 

-g 

The parameter q corresponds to the right end of the Fermi interval [-^;^] on which the ground state's Bethe 
roots condensate. It is fixed by the value of the chemical potential h by demanding that the dressed energy e(/^), 
defined as the unique solution to the below integral equation, vanishes at ±q: 

1 

s{A)- f K{A-u)e{li)^ =£o{A) with £o{A) = A^-h ande(±^) = 0. (2.8) 
J 2n 



' Note that different sets of roots [^(J and [firj lead to different sets of background parameters 

*Tlie uniqueness of solutions follows from the fact that the solution to i2.2l are such that fif^ (0) - /jf^ (0) + 2inf5/L. This allows one to 
show that (a)) is strictly increasing on R + linp/L and maps it onto R. Moreover, one can check that 3 (^4,)) 7^ on C \ (R + 2ot/3/L). 
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We also remind the relation pf - where pf - p{q) is the Fermi momentum. 

In the following, we will focus on the excited states in the N,^ = N + 1 -particle sector only. In order to describe 
the thermodynamic properties of such yS-deformed excited states, it is convenient to introduce the associated shift 
function 



N 



Fif^j (CO) ^f(oj\ {yu^jf ^) =Ly (oj) - ^(41 (cu)] = i- J] e (a; - i,) - i- ^ (a; - - ^ + ^yS . (2.9) 
It can be shown that this counting function admits a thermodynamic limit F/j that solves the linear integral equation 



FA A 



q 



There fip^, resp. jUf,^, are to be understood as the unique solutions to ^{pipj = Pa/L, resp. ^{phj - Ki/L, where 
^ is given by (12.61 ). Note that we have explicitly insisted on the auxiliary dependence of the thermodynamic limit 
of the shift function on the positions of the particles/holes. However, in the following, whenever the value of 
{fipj and {f^hj will be dictated by the context, we will omit it. We also remind that the above shift function 
measures the spacing between the ground state roots Aa and the background parameters ju^ defined by ^jf^j : 

fia -Aa=Fp {Aa) ' [^'{Aa)]-' (l + O (L'^)). 

The integral equation for the thermodynamic limit of the shift function Fp can be solved in terms of the dressed 
phase (p {A, /j.) and dressed charge Z (A) 



4>{A,n)- J K{A-t)4>{t,h)^^ ^9{A-n) 



and 



Z{A) 



J 



K{A-t)Z{t)^ = 1 . (2.10) 
2n 



Namely, 

Fp(A) = F/s\A 



Here, we also remind two very nice relationships that exist between the dressed phase and dressed charge 



Z{A)= l+(p{A,-q)-(p{A,q) and Z~\q) ^ I + (p{-q,q) - (p{q,q) 



(2.11) 



(2.12) 



The first one is easy to obtain and the second one has been obtained in 115811791 . 

The shift function allows one to compute many thermodynamic limits involving the parameters [ni^]- For 
instance, introducing the combination of bare momentum and energy uq{A) = poiA) - tEo{A) /x, one readily sees 
that for a n particle/hole excited state {p.c^ ) at y6 = 



(N+l 



' ~^ °^ ' ■ j^^Q a=\ 



(2.13) 



^0=1 



a=\ 



Above and in the following, u stands for the combination of dressed momenta and energies u(A) = p (A) - ts {A) jx. 
It admits the integral representation 



u (/I) = Mo (/I) 



- Jmq (h) (j) (p.. A) dju . 



(2.14) 
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The function admits a unique zero of first order on R. It is believed tliat tliis property is preserved for u. Clearly, 
in virtue of Rouche's theorem, this holds true for c large enough. We will not purse the discussion of this property 
here as it goes out of the scope of this paper and will use it as a working hypothesis. In other words, we assume 
that given a fixed ratio t/x, there exists a unique Aq such that u'{Ao) = and u"{Ao) < 0. 

We do stress however that this working hypothesis should not be considered as a restriction but a simplification 
of the exposition at most. Indeed, it follows from \u'{A)\ +00 when % (A) ±00 that, for any value of c > 0, 
u' has a finite number of real zeroes. The case when u' has multiple real zeroes of arbitrary order could be treated 
within out method but would make the analysis heavier. 

As a concluding remark to this section, we would like to stress that all functions that have been introduced 
above (the dressed momentum p, the dressed energy e, the dressed charge Z and the dressed phase 0) are holo- 
morphic in the strip 

Us = {zeC : |3(z)|<2(5) (2.15) 

around the real axis. The parameter 6 satisfies c/8 > 6 > and is chosen sufficiently small so that p is injective 
on Us and that one has inf^^u^ [% (Z{A))] > 0. We will tacitly assume such a choice in the following each time 
the strip Ug will be used. 

3 The method and main results 

The zero-temperature one-particle reduced density matrix in finite volume refers to the below ground state expec- 
tation value: 

Pn{x, t) = {ilj (l^lf ) I (D (X, t) 0) I ^ ({^)f ) ) • ({^fllf ^ ■ (3-1) 

The parameters {Aa}^ correspond to the set of Bethe roots parameterizing the ground state of (12.11) . We recall that 
the fields evolve in space and time according to 

^{x, t) = e-'^P+''ff^Ls (D(o, 0) e'-^'P-''^"" , (3.2) 

where Hf^is is the Hamiltonian of the model given in (12.11) and P is the total momentum operator. The action of 
P on the eigenstates of H^ls has been computed in |4|. 

We denote by p(x, t) - limA^j^^+oo Pa?(jc, t) the, presumably existing, thermodynamic limit of pn{x, t). We will 
not develop further on the existence of this limit, and take this as a quite reasonable working hypothesis. 

3.1 Description of the method 

In this article, we carry out several manipulations that lead us to propose a series representation for p{x, t) giving 
a straightforward access to its leading large-distance/long-time asymptotic behavior. 

The starting point of our analysis is the model in finite volume. We will first provide certain re-summation 
formulae for pn{x, t) starting from the form factor expansion of (13.11 ). The latter involves a summation over all 
the excited states {ie over all solutions to (I2.2I )- (I2.3I) at y6 = 0). This sum has a very intricate structure which 
prevents us from analyzing its thermodynamic limit rigorously from the very beginning. We therefore introduce a 
simplifying hypothesis. Namely, denoting the energy of an excited state by and the one of the ground state by 
figs we argue that all contributions issued from excited states such that - £'gs scales with L do not contribute 
to the thermodynamic limit of the form factor expansion of pn{x, t). In the light of these arguments, we are led to 
analyze an effective form factor series pN-efti^, t) and a certain y-deformation pn-&s{x, t \ y) thereof. Our conjecture 
is that PiV;efT(x, t \ y = I) = pN;es{x, t) has the same thermodynamic limit as pn{x, t). 
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We study y i-^ PN-eff i^, ^ I t) by means of its Taylor coefficients at y = 0: 

d'" 

P^iffCx, t) = -^ PN;eS (X, t \ y\^^ . (3.3) 

All rigorous, conjecture-free, results of this paper are relative to these Taylor coefficients. We show that these 
admit a well defined thermodynamic limit p^^{x, t). In addition, we provide two different representations for this 
limit, each being a finite sum of multiple integrals. 

• The first representation is in the spirit of the ones obtained in fl31l63l . It corresponds to some truncation of 
a multidimensional deformation of a Fredholm series for a Fredholm minor. 

• The second representation is structured in such a way that it allows one to read-off straightforwardly the 
first few terms of the asymptotic expansion of p^^{x, t). The various terms appearing in this representation 
are organized in such a way that the identification of those that are negligible {eg exponentially small) in the 
X — > +00 limit is trivial. 

The above two results are derived rigorously without any approximation or additional conjecture. However, 
in order to push the analysis a little further and provide results that would have applications to physics, we need 
to rely on several conjectures. Namely, we assume that 

1 . the series of multiple integrals that arises upon summing up the thermodynamic limits of the Taylor coeffi- 
cients li,t,=oPlff''(-'^' /'w! is convergent; 

2. this sum moreover coincides with the thermodynamic limit of p^-csix, t\y - \) and hence, due to our first 
conjecture, with p {x, t). 

These conjectures allow us to claim that p (x, t) can be represented in terms of a series of multiple integrals. The 
latter series corresponds to a multidimensional deformation of the Natte series expansion for Fredholm minors 
of integrable integral operators |62|. This multidimensional Natte series has all the virtues in respect to the 
computation of the long-time/large-distance asymptotic behavior of p{x, t); it is structured in such a way that one 
readily reads-off from its very form, the sub-leading and the first few leading terms of the asymptoics. 

So as to conclude the description of our method, we would like to stress that the aforementioned conjectures 
of convergence are supported by the fact that they can be proven to hold in the limiting case of a generalized free 
fermion model ||62]| . Unfortunately, the highly coupled nature of the integrands involved in our representations 
does not allow one for any simple check of the convergence properties in the general +co > c > case. 

3.2 Large-distance/long-time asymptotic behavior of the one-particle reduced density matrix 

We have now introduced enough notations so as to be able to present the physically interesting part of our analysis. 

Let X > be large and the ratio xjt is fixed. Let Aq be the associated, presumably unique {cf (12.141) ). saddle- 
point of u{A) - p {A) - tE {A) /x. Assume in addition that Aq t ±q and Aq > -q. Then, under the validity of the 
aforementioned conjectures, the thermodynamic limit of the zero-temperature one-particle reduced density matrix 
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p {x, t) admits the asymptotic expansion 



p'iAo) e-["Wo)-"('?)] \r,'°\ , . 

P^^^^)^ A 777771 77777 77 7 77 — 7 Ul r Uo) + o (1)1 

+ ^(l+o(l)) 

[i{x + VFtW'>'^-'^''-^i [-i{x - VFtW^''^''^i 

+ ^ ^(l+o(l)) + }_/Cf..r— (l+o(l)) (3.4) 

[/ (x + vptW^^-'i^^ [-i {x - VFtW»^'''>*^^ e\t-ez ^ 

The critical exponents governing the algebraic decay in the distance of separation are expressed in terms of the 
thermodynamic Umit F'^'^ of the shift function (at p - 0) associated with an excited state of (12.11 ) having one 
particle at fip and one hole at namely, 

a Z{A) _„ Z{A) >„ Z{A) 



The type of algebraic decay in the explicit terms in (13.41) can be organized in two classes. There is a square 
root power-law decay (te"(/lo) - xp"{Aq))~~- stemming from the saddle-point Aq. All other sources of algebraic 
decay appear in the so-called relativistic combinations x ± Vft and exhibit non-trivial critical exponents driven 
by the shift function of the underlying type of excitation. We recall that ±vf corresponds to the velocity of the 
excitations on the right/left Fermi boundary: vf = e'{q) Ip'iq)- 

Each of the three explicit terms in these asymptotics has its amplitude (|?^''''f , 1^? ^P) given by the 

thermodynamic limit of properly normalized in the length L moduli squared of form factors of the conjugated field 
O ^ More precisely, 

• |?^/°|^ involves the form factor of taken between the A/^-particle ground state and an excited state above 
ih&N + I particle ground state with one particle at Aq and one hole at q. 

I i2 

• \Tq ^1 corresponds to the case when one considers an excited state above the N + 1-particle ground state 
with one particle at -q and one hole at q. 

• corresponds to the case where the form factor average of <1)^ is taken between the N and the N-f1- 
particle ground state. 

The explicit (but rather cumbersome) expressions for the amplitudes together with a more precise definition 
are postponed to appendix IA.3I 

Also, 1] ^ .+oo [ stands for the characteristic function of the interval ] q ; -i-cx3 [. It is there so as to indicate that, to 
the leading order, the contribution stemming from the saddle-point only appears in the space-like regime Aq > q. 
We stress however that hole-type excitations in a vicinity of the saddle-point also contribute in the time-like regime 
where Aq e\-q;q[. This fact follows from the structure of the terms present in the sum over . 

We would now like to discuss the sum over the integers , €~ in (13.41 ). The latter represents the contributions 
to the asymptotics associated to the so-called quicker harmonics. Ineed, every term in this sum oscillates with a 
phase 

ip{.,e- - tu iq) + ru i-q) - {t + r)u{Ao) . (3.6) 
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It is also caracterized by its own critical exponents 



- (1 + r + A+)^ + (A„ - 




(3.7) 



where, 



A± = 



Z{±q) 



- t ^{±q,-q) 



{t + l)(t>{±q,q) + + ncp{±q,Ao) . 



(3.8) 



2 



Our method of analysis only allows us to prove that the only harmonics present in the asymptotics are those 
oscillating with one of the frequencies ipt+/- and that they decay, to the leading order {ie up to o (1) terms), with 
the critical exponent A/'+/-. We are however unable to give an explicit prediction for the amplitudes C{+j-. Note 
that the sum runs over all integers {- subject to the constraint r}{{^ + €~) > 0. The parameter ij depends on the 
regime: = 1 in the space-like regime (Aq > q) and j; = -1 in the time-like regime {\Aq\ < q). Finally, the * in the 
sums indicates that one should not sum up over those integers €' giving rise to the frequencies that are present 
in the first three lines of (13.41) . 

Note that we have organized the large x (with xjt fixed) asymptotic expansion in respect to the various oscil- 
lating phases. Each phase appears with its own exponent driving the power-law decay in x. Our computations 
allowed us to compute the leading {ie up to o (1) corrections) behavior of each harmonic. Note that the o (1) terms 
stemming from one of the harmonics may be dominant even in respect to the leading terms coming from another 
harmonic. 



The oscillating phases and amplitudes appearing in (13.41 ) are reminiscent of the type of excitations that give rise 
to their associated contribution. Each term in (13.41) can be associated with some macroscopic state of the model. 
For instance, the one occurin^ in the first line of (13.41 ) corresponds to a macroscopic state characterized by one 
particle at Aq and one hole at q. There are infinitely many microscopic realizations of such a macroscopic state. 
For instance, any excited state realized as one particle at Aq, one hole at q, 

• n+ particles and holes u, located at q in the thermodynamic limit: u), ,u\ — > qfov a = 1, . . . , 

"a ' "a N,L-^+oa 

• n- particles //p^ and holes located at -q in the thermodynamic limit: yuj,'^ , //|'^ ^ — ^ for a = 1 . 

would give rise to the same (from the point of view of energy £ = s{Aq) , momentum P = p{Ao) - p{q),...) 
macroscopic state. In a joint collaboration with Kitanine, Maillet, Slavnov and Terras we have shown PF3l that 
indeed, in the zero-time case, the contribution of a given macroscopic state to the asymptotics is obtained by 
summing up over all such zero-momentum excitations on each of the Fermi boundaries. Clearly, this picture 
persists in the time-dependent case as well. The only difference being that, in the time-dependent case, the number 
of relevant macroscopic states contributing to the asymptotics is bigger (one has to include the contributions of 
excitations around the saddle-point in addition to the excitations on the Fermi boundary). Moreover, we would like 
to draw the reader's attention to the fact that it is precisely the sum over such zero momentum excitations on the 
Fermi boundary that gives rise, through some intricate microscopic mechanism of summation, to the relativistic 
combinations {x + vpt)"* (in what concerns the left Fermi boundary) and {x - vpt)"' (in what concerns the right 
Fermi boundary) arising in the asymptotics. This mechanism can be considered as yet another manifestation of 
conformal field theory on the level of asymptotics. 

' The o (1) corrections being excluded 



Remarks 
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Our analysis leads us to propose an alternative interpretation of the universality hypothesis. Namely, when 
dealing with asymptotics (large-distance, etc) of correlation functions, one is brought to the analysis of the con- 
tributions of "relevant" saddle-points. As one can expect from the saddle-point type analysis of one-dimensional 
integrals, the leading asymptotics are only depending on the local behavior around the saddle-point of the driving 
term. All other details of the integrand do not matter for fixing the exponent governing the algebraic decay. There- 
fore, it is quite reasonable to expect that models sharing the same types of saddle-points exhibit the same type 
of critical behavior. The universality hypothesis 1*24] stating that models sharing the same symmetry class have 
the same value for their critical exponents can be now re-interpreted as the fact that the symmetries of a model 
uniquely determine the structure of the driving terms in the saddle-points that are relevant for the asymptotics. 
As a consequence, the leading power-law decay stemming from the local analysis around these saddle-points is 
always characterized by the same critical exponents regardless of the fine, model dependent, function content of 
the integrals describing the correlation functions. 

We draw the reader's attention to the fact that the terms appearing in the 2""^ and 3'^'^ lines of (13.41 ) correspond 
solely to excitations on the Fermi boundaries and confirm the CFT/Luttinger liquid-based predictions for the 
long-distance asymptotic^ at t - due to the identifications following from (12.121 ): 

F:(g) + l = ^ , Fli-^) = , F-^(q) = ^-Z(g) , F-^(-q)-l = -^-Z(q) . (3.9) 

However, we do stress that (13.41) clearly shows the need to go beyond the CFT/Luttinger liquid picture so as 
to provide the correct long-time/large-distance asymptotic behavior of the correlation functions in gapless one- 
dimensional quantum Hamiltonians. In particular, our results contain additional terms in respect to the predictions 
obtained in [2|. Our result has a strong structural resemblance with the non-linear Luttinger liquid based predic- 
tions for the edge exponents [28 ] and amplitudes [ 12J arizing in the low momentum k and low energy a> behavior 
of the spectral functioifl 



4 The form factor series 

In this section, we will provide two new representations for the zero-temperature reduced density matrix (13.11 ) 
starting from its form factor issued expansion: 

^ ^ „ n>'"<"-'|(^Kir')K'co.o)k(iA,)f))f 

PN{x,t)= 2j — n ..2 ^ • (4.1) 

The above series runs through all the possible choices of integ ers — l,...,A'^-i- 1 such that £i < • • • < -Ef^^i. 

Below, we shall argue in favor of several reasonable approximations that allow us to reduce the form factor 
series to another, effective one, whose structure is simple enough so as to be able to continue the calculations 
directly on it. 



'Taking the / — > limit of ( I3.4t is sligiitly subtle. The first line produces a contribution proportional to e'T . In the ? — > limit, this 
function approaches, in the sense of distributions, a Dirac 6 (x) function. The presence of this 6 (x) function is expected from the form of 
the commutation relations between the fields. However, in the large-x limit of interest to us, it does not contribute. 

*The latter corresponds to the space and time Fourier transform of <$ (x, t) O"^ (0, 0)>ljo;+oo[ (0 + (O'' (0, 0)0 (x, 0>l]-oo;0[ W 
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4.1 The effective form factors 



It has been shown in ll64ll (slightly different determinant representations for these form factors have already ap- 
peared in ||6T]|73l) that the form factors of the operator taken between the A'^-particle ground state {/la}f and 
any particle-hole type excited state {l^eJi^^ as described in (12.31 ) takes the form 

This representation involves two functionals, the so-called smooth part of the form factor Q^-i and the so-called 
discreet part Dp/. These are functionals of the counting function ^ for the ground state, of the counting function 
^j^^l for the excited state and of the associated shift function F{e„]- 

It has been show in 1.641 . that, in the large L-limit and for any n particle-hole type excited state, with n bounded 
independently of L, these functionals satisfy 

(§naD^) ( ) [F^i.aiuJ] = Qn^ ( l^^^l ) [f 0, ^, ^Fo] • D, ( ^P^"^^ ) [Fo, ^Fo] (l + O (i^)) . (4.3) 

We stress that the functionals appearing on the rhs of the above equation act on 

i) the thermodynamic limit Fq{A) of the shift function at yS = associated to the excited state labeled by the 
set of integers {^if""' (IXTTT) . 

ii) the thermodynamic limit ^ {A) of the counting function (12.61 ). 

iii) the counting function associated with Fq: ^Poi'^) - ^('^) + ^o('^) Z^- 

We do stress that the shift function Fq depends implicitly on the rapidities of the particles {yUp^}" and holes 
entering in the description of the excited state of interest, cf (12.1 II ). We chose not to write this dependence 
explicitly in (14.31 ) as the auxiliary arguments of Fq are undercurrent by those of the functionals Dj^ and 'Qn,\- Given 
any holomorphic function v {A) in a neighborhood of R, the explicit expressions for Da? [v, ^, ^y] (and Qj.^ [v, ^, ^y]) 
involves two sets of parameters {Aa]^ and which are defined as follows 

• iJ.k, k eZis the uniqu^ solution to ^{fik) - k/L, ie the second argument of the functionals; 

• At, k e ^ \ ; N'S is the uniqu^ solution to ^yiAk) = k/L, ie the third argument of the functionals. 

We insist that here and in the following, the parameters fx^ or Ap entering in the explicit expressions for these 
functionals are always to be understood in this way. Also, we remind that the integers £a are obtained from the 
integers {pa}" and {ha}" as explained in (12.31 ). 



• The discreet part 



The functional represents the universal part of the form-factor: 



D 



N 



{haYl 



nti{4sin2[7rFo(/l^)] 

N+1 N 

n 2n^'(M0 n 2K^'^^iAa) 

a=\ a=\ 



dett 



1 



- h 



(4.4) 



'The uniqueness follows from the fact that the dressed momentum p {X) is a biholomorphism on some sufficiently narrow strip Ug 
around the real axis and that p (/i) 6 R => /I 6 R. 

'The uniqueness follows from Rouche's theorem when L is large enough. 
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The large A'^, L behavior of (14.41 ) can be computed explicitly and is given in (IA.2I) - (IA.4I) . However, it is the above 
finite product representation of D^y that is suited for carrying out resummations. 

• The smooth part 

The functional Qnj represents the so-called smooth part of the form factor: 

^^■■y[{hjiy^o,^,^F,l det^^ [^]] det^ [SW) [^^j] ^ {^.j" ) ij [[{v^.^^i^.J ^.^^ - + iec j 



(4.5) 



Above, we have introduced several functions. For any set of generic parameters ({Zo)"; {ya)") £ U'^ x 

. xn \ = \ \ 7 r-7 r- and VN-^ioj) ^ \ — (4.6) 

\ {Jflli / \ \ (ya-Jb- IC) iZa -Zb-lC) \ a)-flb + leC 



Also we have set 



^i^U^^ ^ KiMi^Z^ and S^.f [f^J = 5 - (4.7) 

Finally, for any set of generic parameters ({Zu}"; {ja)"^^) ^ U'^ x the entries of the two determinants in 
the numerator read 



Vkdy]({zayi;{yayr') = 



Y[{zk-ya)Y[{zk-Za + ic) 
.0=1 0=1 a: (xu - Zi) 

n n+1 Q-linv{zk) _ I 

n iZk - Za) n {Zk -ya+ ic) 
ai^k a=l 



n+l n 

_ Ylizk- ya) n izk - Za- ic) 

\Ar^(iz.nMr) = 1^ (4.8) 

n iZk - Za) n iZk -ya- ic) 
a*k a=\ 

Note that the singularities of the associated determinants at Zk - Zj, j i= k are only apparent, cf P31l64ll . 



4.2 Arguments for the effective form factors series 

It is belie veclll that when computing the T - OK form factor expansion of a two-point function (G.S .\O1O2\G.S .) 
on the intermediate excited states (as in (14. Il l), the contribution of those excited states whose energies differ 
macroscopically from the ground state's one {ie by a quantity scaling as some positive power of L) vanishes in 
the L — > -1-00 Umit. This can, for instance, be attributed to an extremely quick oscillation of the phase factors and 
the decay of form factors for states having large excitation momenta and energies. Therefore, we shall assume 
in the following that the only part of the form factor expansion in (14.11) that has a non-vanishing contribution to 

' The computations presented in appendices IB.2l and lB.3l can be seen as a proof of tiiis statement in tlie case of a generalized free 
fermion model. 
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the thermodynamic limit p{x,t) of p^ixj) corresponds to a summation over all those excited states which are 
realized as some finite (in the sense that not scaling with L) number n, n = 0,1, . . . , of particle-hole excitations 
above the (N + l)-particle ground state. Indeed, these are the only excited states that can have a finite {ie not 
scaling with L) energy gap above the ground state in the A'^-particle sector. 

Even when dealing with excited states realized as a finite number n of particle-hole excitations above the 
(A'^ + l)-particle ground state, it is still possible to generate a macroscopically different energy from the one of the 
N-paiticle. ground state if the rapidities of the particles become very large (ie scale with L). This case corresponds, 
among others, to integers pa becoming very large and scaling with L. We will drop the contribution of such excited 
states in the following. 

Limiting the sum over all the excited states in the {N + l)-particle sector to those having the same per-site 
energy that the ground state means that one effectively neglects correcting terms in the lattice size L. It thus seems 
very reasonable to assume that, on the same ground, only the leading large-L asymptotic behavior of the form 
factors will contribute to the thermodynamic limit of p^ (x, t). It is clearly so when focusing on states with a low 
number n of particle/hole excitations. However, in principle, problems could arise when the number n becomes of 
the order of L. Our assumption lead to the following consequences: 

• we discard all summations over the excited states having a too large excitation energy. This means that we 
introduce a "cut-off" in respect to the range of the integers entering in the description of the rapidities of the 
particles. Namely, we assume that the integers pa are restricted to belong to the sejl 

&l;^ = [neZ : -wl < n < wl]\11 ; N + where wl ~ L^*~^ ■ (4.9) 

N+l N 

• The oscillating exponent ^ mq - Z "o i^a) is replaced by its thermodynamic limit as given in (12.131) . 

a=l a=\ 

• We drop the contribution of the O ■ In terms in the large-size behavior of form factors given in (14.31) . 

Note that, within our approximations, the localization of the Bethe roots {yU^lf^^ for an excited state whose 
particles' (resp. holes') rapidities are labeled by the integers {pa}"=i (resp. [ha}"^^) does not depend on the specific 
choice of the excited state one considers. Hence, we effectively recover a description of the excitations that is in 
the spirit of a free fermionic model. 

Our simplifying hypothesis suggest to raise the below conjecture 

Conjecture 4.1 The thermodynamic limit of the reduced density matrix p^ (x, t) coincides with the thermody- 
namic limit of the effective reduced density matrix pN;eS (x, t): 

lim pN{x,t)= lim pNes{x,t) (4.10) 

N,L^+oo N,L—>+oo 

where pN;eS i^, is given by the series 

n=0 Pi<-<P„hi<-<h„ a=l ^ ^ '^ \' «'1/L \ •'^"fll / J 

There Si^ineZ : -wl <n< wl], ^ SlMI; N + \J and S'l^ = l\; N + \J. Also, the * refers to the 
running variable ofpQ on which the two functionals act. 

*Note that we could choose w/, to scale as L'+^, where 6 > is small enough but arbitrary otherwise. We choose e = 1/4 for definiteness. 
cf appendix IB.ll for a better discussion of the origin of such a property. 
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The effective form factor series (14.111 ) possesses several different features in respect to the form factor expansion- 
based series that would appear in a generalized free fermion model (c/ (IB.20I )). Namely, 



• the shift function T^o depends parametrically on the rapidities of the particles and holes entering in the 
description of each excited state one considers, cf (12.1 II) . It is thus summation dependent. 

• Each summand is weighted by the factor Qn-^i that takes into account the more complex structure of the 
scattering and of the scalar products in the interacting model. This introduces a strong coupling between 
the summation variables {pa}'\ and {ha}\. Indeed, the explicit expression for Qt^-\ involves complicated 
functions of the rapidities {^ip„}'l and {yu/i,,}", which, in their turn, depend on the aforementioned integers. 

A separation of variables that would allow one for a resummation of (14.111) is not possible for precisely these 
two reasons. To overcome this problem, we proceed in several steps. First, we introduce a y-deformation of the 
effective form factor series such that pn-cS {x, t \ y)|y=i = Pn;&s i^^ 0- 



„-n D\< -<D,. h.^...^h ^ \riu/ \^ 



n=0 Pi< -<Pn h\<-<h„ a=\ 



'{PaT, 



(4.12) 



For any finite A'^ and L, it is readily checked by using the explicit representations (14.41 ) for and (14.51 ) for 0N;y 
that the y-deformation pN;ef£ (x, t \y) is holomorphic in y belonging to an open neighborhood of the closed unit 
discs Hence, its Taylor series around y = converges up to y = 1. We will then show in theorem fC. II that, given 
any fixed m, the m"^ Taylor coefficient of pN;es{x, f | y) at y = 0: 



d'" 



(4.13) 

7=0 



can be re-summed into a representation where the existence of the thermodynamic limit p^^{x, t) is readily seen. 
This fact is absolutely not-clear on the level of (14.131) as, due to (IA.3I )- (IA.4I ). each individual summand vanishes 
as a complicated power-law in L that depends on the excited state considered. We will then show that one can 
represent the thermodynamic limit p^^{x, f) in another way. This representation is given in terms of a finite sum of 
multiple integrals and corresponds to a truncation of the so-called multidimensional Natte series that we introduce 
below. The latter description of p^"^{x, t) gives a sti^aightforward access to its asymptotic expansion. 

The proof of the existence of the thermodynamic limit and the construction of the truncated multidimensional 
Natte series for p^^{x, f) constitute the rigorous and conjecture free part of our analysis. This is summarized in 
theorem 14.11 

Working on the Taylor coefficients pj^"^j^.(x, t) instead of the full function pN-^six, t \ y) taken at y = 1 has 
the advantage of separating all questions of convergence of the representations we obtain from the question of 
well-definiteness of the various re-summations and deformation procedures that we carry out on p^^^^^^ix, t) (and 

subsequently on p'^^ix, t) once that the thermodynamic limit is taken). Indeed, by taking the y-derivative at 
y = 0, we always end up dealing with a finite number of sums. However, if we had carried out the forthcoming 
re-summation directly on the level of Peff (x, t), we would have ended up with a series of multiple integrals instead 
of a finite sum. The convergence of such a series constitutes a separate question that deserves, in its own right, 
another study. Nonetheless, in the present paper, in order to provide physically interesting results, we will take 
this convergence as a reasonable conjecture in a subsequent part of the paper. 



^The apperent singularity of tlie determinants at e*'"'^°*'''' - 1 = 0, c./ l l4.8t . are candelled by the pre-factors sin^[;TFo(/lt)] present in 
Dn, cf 
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4.3 An operator ordering 

Prior to carrying out the re-summation of the form factor expansion for pj^"gj^(-^, t), we need to discuss a way 
of representing functional translations and generalizations thereof. These objects will allow us to separate the 
variables in the sums occuring in (14.131) . and carry out the various re-summations. A more precise analysis and 
discussion of these constructions is postponed to appendix |Dl In the following, we denote by ^ (W), the ring of 
holomorphic functions in { variables on W c C^. Also, here and in the following f e iff (W), with W non-open 
means that / is a holomorphic function on some open neighborhood of W. Finally, for a set S on which the 
function / is defined we denote \\f\\^ = sup^g^ 1/(^)1- 

Throughout this paper we will deal with various examples {Dn , , • - •) of functionals !F[v] acting on 
holomorphic functions v. The function y will always be defined on some compact subset M of C whereas the 
explicit expression for 'F [v] will only involve the values taken by y on a smaller compacH K c Int (M). In fact, 
all the functionals that we will consider share the regularity property below: 

Definition 4.1 Let M, K be compacts in C such that K c Int (M). Let Wy be a compact in C^-, 4 £ M = {0, 1, . . . }. 

An i ^-parameter family of functionals 'Fl-] (z) depending on a set of auxiliary variables z ^ is said to be 
regular (in respect to the pair (M, K)) if 

i) there exists constants Cj^ > and C > such that for any f,ge0' (M) 

\\f\\K + \\8\\K<CT ^ \\r[f]{-)-r[g]{-)\\„_ <c'\\f-g\\K , (4.14) 

where the ■ indicates that the norm is computed in respect to the set of auxiliary variables z € W^. 

ii) Given any open neighborhood Wy ofO in C^>, for some £y € N, ifv (A, y) € (^M x Wy^ is such that WvH^xWy < 
Cjr, then the function (y,z) ^ 'F[v(*,y)] (z) is holomorphic on Wy x W^. Here, the * indicated the running 
variable A ofv{A,y) on which the functional T [•] (z) acts. 

The constant Cj^ appearing above will be called constant of regularity of the functional. 

This regularity property is at the heart of the aforementioned representation for the functional translation and 
generalizations thereof that we briefly discuss below. However, prior to this discussion we need to define the 
discretization of the boundary of a compact. 

Definition 4.2 Let M be a compact with n holes ( ie C\M has n bounded connected components) and such that dM 
can be realized as a disjoint union of n + 1 smooth Jordan curves ja '. [0 ; 1 ] — > dM, ie dM = U"=i Ta ([0 ; 1 ])■ 
A discretization (of order s) of dM will correspond to a collection of (n + l){s + 2) points tj^a = Jaixj) with 
j = Q, . . . , s + \ and a = \, . . . ,n -\- \ where xq = < xi < • • • < < 1 = x^+i is a partition o/ [0 ; 1 ] of mesh 
21 s: \xj+\ - Xj\ < 21 s. 

4.3.1 Translations 

Suppose that one is given a compact M in C without holes whose boundary is a smooth Jordan curve y : [ ; 1 ] — > 
dM. Let be a compact such that K c Int (M) and T a regular functional (c/ definition 14. II) in respect to (M, K), 
for simplicity, not depending on auxiliary parameters z. 

*Here and in the following, Int (M) stands for the interior of the set M. 
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It is shown in proposition ID . 1 1 that, then, for \y\ small enough one has the identity 



r 



ywJ* jj^'jl j =^lim Ipje^'^J''')-^'^^'') • r[7fs] 



(4.15) 



' lft=0 



The function W„ appearing above is defined in terms of an auxiliary function ij/ {A,fi) that is holomorphic on MxM 



{yaYl 

{Za}1 



Y,^{A,ya)-ilj{A,Za) Whereas /. (a \ {ga]\) ^ fs U) - ^ - A it ' ^"^"^^^ 



Finally, gs (A) is a differential operator in respect to with a = 1, . . . , ^: 



g,iA) = Y,m^)-^.- 



7=1 



(4.17) 



The definitions of and fs involve a set of + 1 discretization points tj of dM. 

The limit in (14.151 ) is uniform in the parameters j„ and Za belonging to M and in \y\ small enough. Actually, the 
magnitude of y depends on the value of the constant of regularity Cf-. If the latter is large enough, one can even 
set y = 1 . The limit in (14.151 ) also holds uniformly in respect to any finite order partial derivative of the auxiliary 
parameters. In particular. 



n I a„P' 



-T 



yWnU 



ly=o 



- lim 



dy" 



kk=o 

\y=0 

(4.18) 




Figure 1: Example of discretized contours. In the Ihs the compact M is located inside of its boundary 'tfout whereas 
the compact K corresponds to K2q as defined in (14.271) . In this case M has no holes. In the rhs the compact M is 
delimited by the two Jordan curves 'rfj„ and ^ouf depicted in solid lines. The associated compact K (of definition 
14.11 ) corresponds to the loop 'tf (Ka) depicted by dotted lines. The compact M depicted in the rhs has one hole. 
This hole contains a compact Ka inside. 



We refer to appendix |D] for a proof of the above statement. Here, we would like to describe in words how 
formula (14.151 ) works. By properly tuning the value of y and invoking the regularity property of the functional 
f^iyfs] one gets that, for any s, {^aVi ^ J^iyfs] is holomorphic in a sufficiently large neighborhood of € C*. 
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This allows one to act with the translation operators 0^=1 e^'*^*'''^"*''^^*^ Their action replaces each variable q 
occurring in by the combination ^^1=1 [^(ta,yb) - il^ita,Zb)]- Taking the limit s +00 changes the sum over 
ta occurring in /, into a contour integral over 'Wout, cf Ihs of Fig. [T] Due to the presence of a pole ait - A, this 
contour integral exactly reproduces the function Wn that appears in the rhs of (14.151) . 

Note that such a realization of the functional translation can also be build in the case of compacts M having 
several holes as depicted in the rhs of Fig. [U Also, there is no problem to consider regular functional [•] (z) 
that depend on auxiliary sets of parameters z- 



4.3.2 Generalization of translations 

In the course of our analysis, in addition to dealing with functional translations as defined above, we will also 
have to manipulate more involved expressions involving series of partial derivatives. Namely, assume that one is 
given a regular functional 'F[f,g\ of two arguments / and g. Then, the expression : d'^f [yfs,'gs]\y=Q '• is to be 
understood as the left substitution of the various d^^ derivatives symbols stemming fromg^ . 
More precisely, let be the below holomorphic function of ai, . . . , 

.V 

g,U) = Y,^itj,A)aj . (4.19) 

The regularity of the functional 'F ensures that the function {op} 1-^ dyF ['yfs,'gs] is holomorphic in ai, . . . , 
small enough. As a consequence, the below multi-dimensional series is convergent for aj small enough: 



We stress that as fs (14.161 ) is a holomorphic function of ^1, ... , ^.y, the functional of coefficients of the above 
series give rise to a family of holomorphic functions in the variables gi,. . .,gs- This analyticity follows, again, 
from the regularity of the functional F [f, g] and the smallness of I7I. 

The : • : ordering constitutes in substituting aj "-^ d^., j = \,. . . ,sm such a way that all differential operators 
appear to the left. That is to say, 

:^-rb/„?.U:. ^]^\^—^.X\\^-—y-^^^^^^^ . ,4.21) 

Although there where no convergence issues on the level of expansion (14.201 ). these can a priori arise on the level 
of the rhs in (14.211) . Clearly, convergence depends on the precise form of the functional T, and should thus be 
studied on a case-by-case basis. However, in the case of interest to us, this will not be a problem due to the quite 
specific class of functional that we will deal with. 
At this point, two observations are in order. 

• (14.211 ) bears a strong resemblance with an 5'-dimensional Lagrange series. 

• The functional (of /v) coefficients appearing in the rhs of (14.211 ) are completely determined by the func- 
tional 'F[yfs,^s\ whose expression only involves standard {ie non-operator valued) functions. Should this 
functional have two (or more) equivalent representations, then any one of them can be used as a starting 
point for computing the coefficients in (14.201 ) and then carrying out the substitution (I4.21I) . 



19 



Actually, for the class of functionals that we focus on, no convergence issues arise. Indeed, in all of the cases, 
the m}-^ y-derivative at y = of the : • : ordered functionals of interest appears as a finite linear combinations (or 
integrals thereof) of expressions of the type 



== ^ n ■ n ^"""'^■"'^ ■ ^ t^^^] = "^^^'^ a« € E 1 ; A^l and e„, € {±1} . (4.22) 

U=l b=\ j|y=0 



Above ya are some auxiliary and generic parameters whereas Aa„ are implicit functions of y and ^i, . . . , gs- For 
L-large enough, Aa„ is the unique solution to the equation ^yf^^a^) - 0!a/L. 

The prescription that we have agreed upon implies that one should first substitute "J^^s "-^ 'gs as defined in (14.191 ). 
Then, one computes the m"^ y-derivative at y = of (14.221 ). this in the presence of non-operator valued functions 
'gs. In the process, one has to differentiate in respect to y the functional 9^ [yfs] and the arguments of i^aj- 
Using that ia^i^^o = l^a„, one arrives to 



^ ^ n ^^'''^'"'^ n ■ ^ = n n ^"'''^'"^ z n ■ ^{".i ^^^^ ■ ^^^-^^^ 

^ U=l b=i l\y=o b=l b=l ni,...,n,=0 j=l 



The sum is truncated at most at iij = m, j - I, ... ,m due to taking the y-derivative at y = 0. It is readily 
verified that the {ny)-dependent coefficients Cj„ j [fs] are regular functionals of fs with sufficiently large constants 
of regularity. It remains to impose the operator substitution on the level of (14.231) aj d^. with all differential 
operators (9^^ , A: = 1 , . . . , 5 appearing to the left. It is clearly not a problem to impose such an operator order 
on the level of the polynomial part of the above expression. Indeed, the regularity of the functionals Cj„.j [fs] 
implies that these are holomorphic in ^i,...,^^ belonging to an open neighborhood A^o of e C*. Hence, 
n^=i ^ft* ■ C{„^ j [/v]|q=o is well-defined for any set of integers {mu). In fact, in all the cases of interest for us, the 
neighborhood A^o is always large enough so as to make the Taylor series issued from the products of translation 
operators YYa^^ Q<^ags{iiaa) gf^igiCv*) convergent. Their action can then be incorporated by a re-definition of 



leading to 

m s 
«i,...As=0 ;=1 '^^j 

with 



Z n{57}-|'v}[/4«^o (4-24) 



fs{X}=fs{X) + y j y {tt,n,,) + y (?,, y,) I . 

In this way, one obtains a (truncated to a finite number of terms) s-dimensional Lagrange series. The procedure 
for dealing with such series and taking their s +oo limits is described in proposition ID.2I In the following, all 
operator valued expressions ordered by : • : should be understood in this way. 



4.4 Resummation of the finite- volume Taylor coefficients 

In order to carry out the re-summation of the effective form factor expansion with the help of functional translations 
and generalizations thereof, we need to regularize the expression for the functional 0N;y with the help of an 
additional parameter /?. This regularization will allow us to represent it as a regular functional that, moreover, has 
a form suitable for carrying out the intermediate calculations. 
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The parameter 

It is easy to see that 

(5;vfc)(j^;;{I)[ri^o ;^ ■,^yF,]=j^^\D^{^p^Jj)^[yFp ; ; v.]} (4-25) 

We now introduce a prescription for taking the j3 ^ limit. When considered as a separate object from D^, 
the functional Q^-y may exhibit singularities should it happen that y'^ {e^"'y^p'-'^j^ _ ij - q, c/ (I4.5I )- (I4.8I ). For \y\ 
small enough, as it will always be the case for us, such potential zeroes correspond to the existence of solutions to 
F^{Aj) = 0. For /3 € U/j^ with 

Ufi, = {zeC : 10 51 (jSo) >^{z)>% (jSo) and |3 (z)| < 3 (j3o)] (4.26) 
51 (;8o) > large enough and 3 (fio) > small enough, there are no solutions of 



1 = for ojeUs, this uniformly in < « < m and (j3, {ya}"^, {za}") e U^, x U'^ x [/^ 



It is clear that the optimal value of ySo preventing the existence of such solutions depends on the width S of the 
strip Us and on the integer m. 

Hence, our strategy is as follows. We will always start our computations on a representation that is holomor- 
phic in the half-plane 2lyS > 0, as for instance (I4.12I )- (I4.13I ). In the intermediate calculations whose purpose is to 
allow one to relate the initial representation to another one, we will assume that y6 € t/^^. This will allow us to 
avoid the problem of the aforementioned poles and represent Qn;-/ in terms of a regular functional that is moreover 
fit for carrying out the intermediate calculations. Then, once that we obtain the final expression, we will check 
that this new representation is in fact holomorphic in the half -plane > and has thus a unique extension 

from f/^Q up to y6 = 0. As the same property holds for the initial representation, both will be equal atp = 0. 

Having agreed on such a prescription for dealing with the y6-regularization and treating the y6 — > limit, the 
effective form factor expansion-based representation for p^"^jj(x, t) (14.121 ) can be simplified with the use of the two 
properties below. 

The functional QN;y 

Given A e R^, we define the compact Ka contained in Us- 

KA^{zeC : |3z| < 6 , |5lz| < a) , (4.27) 

and denote the open disk of radius r by 2)o,r = {z e C : |z| < r). 

As follows from lemma IA!21 given A > and large and m e N* fixed, there exists 

• a complex number /?o with a sufficiently large real part and an imaginary part small enough 

• a positive number yo > small enough 

• a regular functional '^^^ 
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such that, uniformly in < « < m, {fipj", ilUhJ") e 2)o,7o ^ ^/So x ic!^ x one has 



HI* 



with // A 



^ 1 1 



The * in the argument of "^^l appearing above indicates the running variable of H on which this functional acts. 

The explicit expression for the functional ^^^^ is given in lemma lATTl The main advantage of such a representation 
is that all the dependence on the auxiliary parameters is now solely contained in the function H given in (14.281) . 
The constant yo is such that 

ji J < i uniformly in (y,/?, {ya}\, {zj\) e £>o,yo ^Up,xKlxKl scuA <n <m . (4.29) 

The functional "^^l is regular in respect to the to the pair {M<^^ , (Ka)) where ^ (Ka) in a loop in Us around 
Ka as depicted in the rhs of Fig.[T]and M<^^ corresponds to the compact with one hole that is delimited by and 
'tfout- This hole contains Ka- Finally, the parameters /Sq eC and yo > are such that the constant of regularity C^^ 
of 1/^] satisfies to the estimates 



nd{dM^^,'W{KA)) 
|5M<yJ + 2nd{dM<f^,'^{KA)) 



>A, 



(4.30) 



where \dMc^^ \ stands for the length of the boundary 5M^^ and d{dMc^^,'^{KA)) > stands for the distance of 
"if (Ka) tod M.^^. 

Similarly to the discussion carried out in section 14.3.21 and according to proposition ID.ll one has that, uni- 
formly in n, p € {0, . . . , m], and Zj, yj, j = I, ■ ■ ■ ,m belonging to Ka'. 



dyP 



H\* 



lr=o 



^ 7=1 



7=0 



(4.31) 



The compact M^^ has one hole. Hence, as discussed in section 1531 one has to consider two sets of discretiza- 
tion points fi p, p - \, . . . ,r + \ fox %„ and t2^p, p = 1, . . . , r -i- 1 for 'ifout- The function m,- appearing in (14.311 ) is 
a linear polynomial in the variables rja^p with a = 1, 2 and p = 1, . . . , r: 



I {ria,p]) - y r 

p=\ 2in\ti^p-A) 



m,p + y ] r 

p=i 2m[t2,p-Aj 



(4.32) 



Finally, g2,r W is a differential operator in respect to rja^p with a = 1, 2 and p = I, . . . ,r: 

d ^ I d 



82. 



p=l ^^'P 



tip - A drji^p ^ t2,p - A di]2,p 



(4.33) 



The functional Dfj 

One can draw a small loop 'tfo^t around K2q in Ug as depicted in the Ihs of Fig. [T] Let be the compact without 
holes whose boundary is delimited by %,ut- Then, given L large enough, the functional D^, as defined by (I4.4I ). is 
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a regular functional (in respect to the pair (Mg, K2q)) of yFp with j3 e U^^ and lyl < yo- The parameters ySo and 70 
are as defined previously. This regularity is readily seen by writing down the integral representation: 



\,...,N (4.34) 



which holds provided that L is large enough (indeed then all /l/s are located in a very small vicinity of the interval 
[-q;q\). Therefore, according to the results developed in appendix IdI and outlined in section 1431 one has that, 
uniformly in y6 € f/gg and < p,n < m 



= lim 

»+oo 



lr=o 

n 

e?l,i(A'pfl)-?l.i(Wa 



a=\ 



•).|^{5«(«)[r.:f:.„,] 



7=0 



(4.35) 



by 



The function A 1-^ v., {A) appearing above is holomorphic in some open neighborhood of K2q in and given 



(a I ^ V, (i) = (^y? - 1/2) Z (i) - (^, q) + y (^^^l-J^l . SL 

^ ' ^—^ ti-A 2m 



(4.36) 



The parameters tj, j = I, . . ., s correspond to a discretisation (cf definition 14.21 ) of the loop ^j^t around in 
Us that has been depicted in the Ihs of Fig. [T] ^-y are some sufficiently small complex numbers and §"1 ^ (A) is a 
differential operator in respect to ga- 



d 



(4.37) 



We remind that the pai^ameters Aa appearing in the second line of (14.351 ) through the expression (14.41 ) for D]^, 
are the uniqu^ solutions to ^yy^ {Aa) - a/L. As such, the /i„'s become holomorphic functions of {^aVi when these 
belong to a sufficiently small neighborhood of the origin in C*. 



Representation for the Taylor coefficients 

To implement the simplifications induced by the functional translations on the level of p^"^jj(x, t), we first observe 
that all of the rapidities fip^ and occurring in the course of summation in (14.121 ) belong to the interval [ -A^ ; Bl ] 
with L^{-Al) = -wl - 1/2 and L^{Bl) = wl + 1/2 (Al > Bi). Hence, a fortiori, they belong to the compact 
K2Al- We can thus represent the smooth part functional as '^^2Al interested solely in the y-derivative 

of (14.121 ) at y = 0. As DN{{paY{,{ha]\) « y^^""'^ and VujA has no singularities around y = Q, all terms 

issuing from n particle/hole excitations with n > m will not contribute to the value of the derivative. Hence, 

^Here, as previously, the uniqueness follows from Rouche's theorem. By writing down an integral representation for ^^l^, one readily 
convinces oneself that, for 7 small enough and given any fixed s, is holomorphic in {^aVi- It is also holomorphic in 7 belonging to some 
open neighborhood of y = 0. 
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we can truncate the sum over n in (14.121 ) ntn - m. Once that the sum is truncated, we represent the functional 



\Dm ■ Sf,'^^ ,1, „ with the help of identities (1435] ) and (|43B- This leads to 



(m) ^ .\ 



lim lim lim 

yS— >0 .V— ^+00 r— >+oo 



SSI: 

n=0 Pt<-<Pn hi<-<h„ a=\ ^-^t^PrJ 



dy" 



I 7=0 



We have set 



Ef (A) 



_ -ixu{A)-g{A) 



with g (A) = gi,, (A) + g2,r (A) 



(4.38) 



(4.39) 



Above, in order to lighten the notation we have not written explicitly the dependence of v^, rur on the auxiliary 
parameters gp, ria,p nor the one of E- {A) on the discretization indices r and s. However, we have kept the hat so as 
to insist on the operator nature of E . We do insist that (14.381 ) has to be understood as it was discussed in section 



Starting from representation (14.381 ). pJ^^'flC-^. can be related with the -derivative of the form factor like 

representation of the functionajl] |^'yVi,£'^j given in (IB. 201 ). Namely, for such an identification to hold, one 
has to extend the upper bound in the summation over n from m xvp to N + \. This does not alter the result as it 
corresponds to adding up a finite amount of terms that are zero due to the presence of y-derivatives. Then, one 
should use the identity 



dy" 



7=0 

?l,=0=Va.p 



?l,=0=Va,i, 



Just as it is the case for the parameters Aj appearing in the expression for D^, the ones appearing in the pre-factors 
of the rhs in (14.401 ) are the unique solutions to ^yv.,(^,s) - s/L. (14.401 ) is an expression of the type (I4.22I ). and to 
deal correctly with it one should implement a : • : prescription for the way the diff'erential operators 5^-^ or d,^^ ^ 
should be substituted in the rhs of (14.401) . 

With the help of identity (14.401) . one is able to force the appearance of the product of function whose 
presence is necessary for identifying the sum over the particle-hole type labeling of integers in (14.381) with the 
functional d'!^XN \ yv^ E^ given in (IB. 201) . This leads to the below representation: 



o^""^ (x t) 



lim lim lim : - — 

/3-*0 i— >+oo r— »+oo oy"' 



(4.41) 



4.5 Taking the thermodynamic limit 

It is shown in appendixO theorem fC. II that p^™^„.(x, t) admits a well defined thermodynamic limit that we denote 



p'^ix, t). This limit is given in terms of a multidimensional analogue of a (truncated) Fredholm series. This series 



^The latter is a functional of v, and g as discussed in subsection l4.3.2l 
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is close in spirit to the type of series that have appeared in B31l63]| . It is also shown in that appendix (proposition 
IC.ll) that it is allowed to exchange 

• the thermodynamic limit N,L ^ +00, N/L D 
with 

• the dy differentiation along with its associated operator substitution, 

• the computation of the translation generated by^2,f-> 

• the computation of the ^-dimensional Lagrange series associated with'gi s, 

• the computation of the r +00 and s — > +00 limits, 

• the analytic continuation in yS from £/^q up to/3 = 0. 

The result of such an exchange of symbols is that p^^-(x, t) admits the representation 

1 

pf^{x,t)= lim lim lim lim : — J £_2(^) . e ' X (.Jry,,^^] c^f^l [^^] U .(4.42) 

This formula deserves a few comments. In the case of complex valued functions £"_, the functional X,Ay,\yVs, 
appearing in (14.421) corresponds to a Fredholm minor (IB. 341) of an integrable integral operator I -\- V acting on 
([-<?;<?])■ The kernel V of this operator is given by (IB. 351) . 




Figure 2: The contour consists of the solid line. The contour 'io^ corresponds to the union of the solid and 
dotted lines. The localization of the saddle-point Aq corresponds to the space-like regime. Both contours lie in 

Usii. 



The subscript '^^"'^ in X,.M[yVs,E^^ refers to an auxiliary compact contour entering in the definition of the 
kernel V. The parameter w delimiting the size of this contour plays the role of a regularization. The limit of an 
unbounded contour "lo^^^ can only be taken after r and s are sent to infinity and the analytic continuation up to 
y6 = is carried out. Finally, in (14.421) also appears the functional It 

can be thought of as the thermodynamic 

limit of the functional '^2w- precise expression and properties are discussed in lemma IATTI 

We also would hke to stress that the parameter /?o defining the region Up^^ from which one should carry out 
the analytic continuation up tOyS = depends on 2w as stated in lemma IATT] This dependence is chosen in such a 
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way that the constant of regularity C^^ for the functional ^^2w l^g^ enough so as to make licit all the necessary 
manipulations with the translation operators and generalizations thereof. 

We stress that formula (14.421 ) constitutes the most important result of appendix O Indeed, it provides one 
with a convenient representation for the thermodynamic limit p^^{x, t). The latter constitutes the first step towards 
extracting the large-distance x and long-time t asymptotic behavior of p^"^{x, t). The proof of such a representation 
for the thermodynamic limit is however quite technical and lengthy. It can definitely be skipped on a first reading. 
Moreover, should one be solely interested in a "short path" to extracting the asymptotics, we stress that formula 
(14.421 ) can be readily obtained without the use of any complicated and long computations. It is enough to take the 
thermodynamic limit formally on the level of formula (14.411 ). Such a formal manipulation leads straightforwardly 
to the representation (14.421 ). 



4.6 The multidimensional Natte series and asymptotics 

Theorem 4.1 The thermodynamic lin 
mensional Natte series representation 



Theorem 4.1 The thermodynamic limit of the Taylor coefficients p'^{x,t) admits the below truncated multidi- 



_ + io^M^M' (X + tvF^y^'^-^'>-'f {x - tvp + /o+)!^^«(^)+'1' (x + tvF)l^'t(-4 

V f Hl}f{{u{z,)}-{zt})[yFl:nyFt;p] (q / {z,} \ ^1 
There, we have introduced the notations 

{z+l = [zt,te : e, = 1} , [z-] = [zt , t e : = -1} , ||z+|| = ^zt , t e : q = 1} . (4.44) 
Fg, Fq \ F^ have been defined in (13.51 ) and, in general, we agree upon 

Ft U) ^f[a ) = - 2] zt) + J] zt) . (4.45) 

6, = 1 £, = -! 

The function ^1,*?^ is related to the thermodynamic limit of the smooth part of the form factor Its expression can 
be found in (IA.8I) . The functionals S, Jl± and are given by 

«r 1 [V] i-g)r-'' G\l + y{q))G\l - v(-^))e-f(-^'^)-^-^)) U/^^^'^^^^^^ 

!t>\v p\ = e ' (4 4o ) 

M (q)]''"^ [2qp'iq)f '^^ [2qp'i-q)f -'^^ (Iny^l^-^'^-^^ 



where G is the Barnes double Gamma function, 
-2qx-^[v]iq) ^1 1 + v(^) \ _ 



-2qx-^[v]{q) ^1 \+v{q)\ 1 / f v{A)-v(m) , 
^+ = ,.7vr.wi r^ I I .-9,w.^ . ' M (^) = exp ^ - J^^ ^—^1^ \ > (4-47) 
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and 

^-[^'^]^^^^( v(V) e-^'-t'-i :^o[v]=e-'..-MUo)(^) - (4.48) 

The second sum appearing in the last line of (14.431) runs through all the elements k belonging to 

TCn - \ k = (ki, . . . , kn+i) : k„+i e M* and ka eN ,a - I, . . . ,n such that ^ aka + k„+i - n \ . (4.49) 
Once that an element ofKn has been fixed, one defines the associated set of triplets J^^y 

/j^-j - { (^1, t2,t^) ,tiel\; n+n,t2el\; kt J, hell; t^- n6,,,„^i j} . (4.50) 
The third sum runs through all the elements {et}tejyj^^ belonging to the set 

B„(Jc) = I {e<}<s/|-, : et e {±1, 0) e J^j^^ with J^et^O for - 1, . . . , « and ^ = 1 i . 

^ t3 = l p=l ^ 

In other words, Bn{k) consists of n-uples of parameters €t labeled by triplets t = {ti,t2,t^) belonging to J^j^y 
Each element of such an n-uple takes its values in {±1, 0). In addition, the components of this n-uple are subject to 
summation constraints. These hold for any value oft[ or t2 and are different whether one deals with ti = \, . . . ,n 
or with tl - n + I. 

The integral appearing in the summand occurring in the third line of (14.431) is n-fold. The contours of 
integration '^^'^^ depend on the choices of elements in &n{k ) and are realized as n-fold Cartesian products of one- 
dimensional compact curves that correspond to various deformations of the base curve 'rf^'' depicted in Fig. |2] 
In the w +00 limit, these curves go to analogous deformations of the base curve '^^\ All these contours lie in 
Usil 

The integrand Ii^n\x^ ({w(z<)l ; [Zt)) M is a regular functional of v, that is simultaneously a function of u(zt) 
and Zt with t running through the set J^j^y This functional depends on the choice of an element {e<}<g/ . from £>n{k ) 
and on x. It appears originally as a building block of the Natte series (cf appendix \B. 5\ for more details). 

We stress that all summands involving the functional are well defined a.t/3 - 0. The potential singularities 
present in are canceled by the zeroes of the pre-factors. 

Proof — 

As a starting point for the proof, we need to introduce the below set of functions depending on the auxiliary 
parameters ap, bi p and Z?2,p- As it has been discussed in section |4.3.2[ these functions will allow us to compute 
the (functional) coefficients necessary for carrying out the operator substitution. We set 

where 

.V ' b ' b 

lu (^) = -y 4>{tp, X)ap and 'g2,r (^) - V + V - (4.52) 
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It is readily checked with the help of lemma IATT] and proposition IB . 31 that for y small enough T given below is a 
regular functional of v,,, vjy and g: 

- /['•^«'(^)+I'(^)]rv,(i)di _ 
T \yy,,g, Tnr\y) = E^{q) ■ e X^.v,[yv„ E^] ^^^Jrur] ■ (4.53) 

In particular [yvs,'g, Wr](y) is holomorphic in y, at least for y small enough. In order to implement the operator 
substitution, we have to compute the Taylor coefficients of the series expansion of [yvs,^, tfrJCy) into powers 
of bi^p, b2,p with p = 1, . . . , r and ap with p = I, . . ., s. These Taylor coefficients are solely determined by the 
functional !F [yv^, g, rur](y) depending on the classical function g ( 14.191 ). Therefore, one can use any equivalent 
representation for 'F [yvs,'g, Wr](T) as a starting point for computing the various partial derivatives in respect to 
bj^p or ap. In other words, one can use any equivalent series representatioi|l]for the Fredholm minor X,^m [yv^, E^]. 
Clearly, different series representations for the Fredholm minor will lead to different type of expressions for the 
Taylor coefficients. However, in virtue of the uniqueness of the Taylor coefficients, their values coincide. As 
shown in [62l, the Fredholm minor we're interested in admits the so-called Natte series representation. The latter 
series of multiple integrals is built in such a way that it gives a quasi-direct access to the asymptotic behavior of 
X^w)[yv, It is thus clear that this is THE series representation that is fit for providing the large-distance/long- 
time asymptotic expansion of the two-point function. We will thus take this series representation as a starting 
point for our calculations. 

The first remarkable consequence of the use of the Natte series is that the exponential pre-factor in front of 
Xc^(n)[yvs, E]:] in (14.531 ) exactly compensates the one appearing in the Natte series (IB.48I ). Once that these pre- 

factors are simplified, one should take the y-derivative of the remaining part of the Natte series representation 
(IB.48I) for X<^(H) [yvs, E^]'^^^^ i^r]- One of the consequences of taking the m"^-y derivative is that the Natte series 

given in (IB.48I) becomes truncated atn = m due to the property //) of the functions i^f appendix IB.5I) : 



^ ,iAui-,)-ui,^,ii-,)-ii,)^ f (g^-)[yv.;^ + /o^] ..(/.) 1 ^ aF_ f (g:^.)[yv.;. + /o^] .^(,) 



(4.54) 



It follows from lemma lA.ll representation (IB. 501 ) and the explicit formulae for the functionals S, and Jl± 
(I4.46I )- (I4.48I ) that the functionals occurring in (14.541 ) are all regular (cf definition 14. II ). Moreover, as follows from 
the previous discussion relative to the procedure of taking the yS — > limit, at this stage of the calculations, 
21 (yS) > is large enough so that the constant of regularity C%^^, of the functionals '^^2w sufficiently large to be 
able to apply proposition ID. 1 1 and corollary ID. II (due to the estimates (14.301) for Cc^jw' constant yo occurring in 
(ID.4I ) is greater then 1 for w large enough, which is the limit of interest) to this functional. Proposition ID. 1 1 and 
corollary ID. II are also directly applicable to all functionals of yv^ in as much as, at the end of the day, one sets 
y = 0. 



^One natural representation that can be used as a starting point for taking the derivatives is the Fredholm series-like representation for 
X^(„.)[7Vj, E^]. In fact, it is this series representation that has been used for the computations carried out in theroem rC.il 
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Clearly, there is no problem to implement the substitution Up i-^ and bi^p i-^ d,^. on the level of (14.541 ) in 
such a way that all the partial derivative operators appear to the left of all r]i p and gp dependent functions. The first 
two lines in (14.541 ) will give rise to translation operators. In the case of the ultimate line in (I4.54I ). this operator 
substitution will produce expressions of the type 



z2]n|i5^[nn|^S7[^nK <k,i)i"' n n 

i„>0«„„>Op=l V^P- ^^P } p=l a=l ['''^"'P'- Or]a,p ) J „=i „=i 



1 d""" 



np>0 «a,p>0 p=l 



dy" 




p=l a=\ 



d"z, 

7=0 {linf 

|?p=0=;7a.p 



(4.55) 



Where Q.p and ^ take the form 



^p{{Zt}) = <p(tp,q)- Ye,(t>(tp,zt) and ^{[z,}) = - — ^— + V, ^t—^ 



Zt 



(4.56) 



One can compute the r, 5 — > +oo limit of such series of integrals by applying corollar#l lD. II and observing 
that ^i;'^ is a Cartesian product of a finite number of compact one dimensional curves that are contained in Us/2- 
In fact, the result of this corollary allows one to cany out the operator substitution in (14.541 ) directly under the 
integration sign. In other words, one is allowed to replace g^i ^ "-^ §"1.5 and ^2,r "-^ 'g2,r directly on the level of 
(14.54 1 ), this without pulling out the partial gp or t/q ^ derivatives out of the integrals. Hence, one is brought to 
computing the action of translation operators. The latter can be estimated by applying proposition ID. II Again, 
there is no problem to apply this proposition either because we compute the y-derivative at y = (so that y 
can be as small as desired in the case of functional of yv,) or because the constant of regularity is large enough 
for '^^2w- follows from this proposition, one can permute the partial y-derivative symbols at y = with the 
action of the finite s and r translation operators. It then remains to take the r +00 and the s — > +00 limits. As in 
each case the convergence is uniform, the limit can be taken directly under the finite sum, compact integrals and 
partial y-derivatives symbols. 

Then, in order to compute the effect of the s — > +cx3 limit we apply the identity (14.151 ) (also cf appendix ID. 31 ): 



lim 



m ]-[' 
+00 1 i 

a=l 



.gl,j(Zo)-gl,iO'a) 



^[yy,] = ^[yF/i] with F/siA) = FJA 



{Za}1 

{yan 



(4.57) 



valid for any regular functional |y| small enough and Za, Ja all lying in IJ^- Here, we would like to remind that 
Fp, appearing above corresponds to the thermodynamic hmit of the j6-deformed shift function, cf (12.111 ). Similarly, 



lim rfe^' ' 

r— >+oo i i 
a=\ 



2.r(Za)-g2,r()'a) , 



■y;2w y',2w 



Hi* 



{ZaYl 

{yan 



with HI A 



izan 

{ya}1 



" 1 



1 



A-Za A-Ja 

(4.58) 



All this for {{ZaYl {ja)") £ K^'^. Then, by applying lemma lATTl backwards. we get 



y;2w 



H\* 



{zj'l 

{yaYl 



(4.59) 



This corollary can be applied to 5f'5„, precisely because its constant of regularity is large enough. 
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The function Q^l, has been defined in (IA.8I) . 



Therefore, we obtain 



grx[n(-g)-M(5)] , _ 



x//«?'\{«fo)l;|z.l)[r^-]C.,.( {J}u{,} (4.60) 

Here F^! i^) ^ Pt (^) + C-^) and f has been defined in 

Once that the functional translations have been computed, one should carry out the analytic continuation of 
the expression in brackets from j6 e f/go up to /J = and then send w to +oo. For this, we recall that the functions 
H^n\x^ admit the below decomposition (c/ (IB.50I )): 

^^?'\{«fo)|;|za)[Al-^l?"({r^t!fo)},{«fc)),to)) W (k[fzw,))^'' f] (e-^-rfj^fe) _ ij' . (4.61) 

It follows from the way Z/,^,-^''^ depends on the set of its v-type arguments (14.611 ) and from the expression for the 
functional 3 [v, u\ KA^ and g^}y ^KM that aU of the expressions one deals with contain the combination 

g2 (i - yFli i-g)) G^{l+ jFt iq)) n (e-2-^^- ^'^ - l)' det.^,^ [/ + YViyFt]] det^^^,, [/ + YViyFt]] 

In virtue of proposition lA.li the function appearing above is holomorphic in {fi,y,{z+},{z-}) € {!R(,6) > 0) x 

The function H^^^^^ {{yFl*_ {zt)},{u{Zt)},{zt}) is analytic in (y,y6) e "Dq^^ x Up^ (here yo is chosen so that 
\yF\l {zt)\ < 112 uniformly in the variables zt, t € J^^^ belonging to Us), and integrable in respect to the {zAtej^,,^^ 

variables. The remaining part of has also the same properties. As the integrals are compactly supported 
it follows that the whole expression appearing inside of the "big" brackets in (14.601 ) is holomorphic in {y,P) € 
Do.yo xU/jg. As a consequence, the -derivative at y = can be continued up to yS = 0. To get the value of the 
analytic continuation at this point it is in fact enough to set y6 = in (14.601 ). 

The last step consists in taking the limit w +oo. This operation will result in an extension of the integration 
contours from bounded ones ^^^'^ to ones going to infinity "if^^K Hence, one needs to check that the resulting 
integrals will be convergent. Note that the function F^'^ (zt) are bounded whenever zt or any of the variables 
belonging to the set {z.+} or goes to infinity. Also, the function ^J^j is bounded at infinity by a polyno- 
mial in Zt of degree n, this uniformly in respect to y-derivatives of order 0, . . . , m. Therefore, as the functions 
^n;v'^ {i^z- iZt)]AuiZt)},{zt}) go to zero exponentially fast in all directions where 'ifj;^-' goes to oo, the integrals 
over 'tf^^'' are indeed convergent. 
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4.7 Some more conjectures leading to the dominant asymptotics of p{x, t) 
Under the assumption that 

1. the Taylor series Y^^q y"^P^^ix> /"j! is convergent up to y = 1, 

2. its sum gives p{x, t), 

3. the multidimensional Natte series given below is convergent. 

We get that p{x,t) is obtained from (14.431 ) by removing the y-derivative symbol and setting y = 1. It then 
remains to identify the coefficients in the first two lines with the properly normalized thermodynamic limit of form 
factors of the field as given in (IA.46I ). (IA.47I ) and (IA.48I ). One then obtains the below series of multiple integral 
representation for the thermodynamic limit of the one-particle reduced density matrix: 



p{x,t)= J X — — — -jl-i ^ooK/to) 

^ \^ a \ K 01 



It follows from the above representation and from conj ecture IB . 1 1 that 
Corollary 4.1 The reduced density matrix admits the asymptotic expansion as given in subsection 



Proof — 

The proof is immediate as far as the multidimensional Natte series defining p (x, t) is convergent. Indeed, 
then, the fine structure of the functions //jj'^''^ given in (IB. 511 ) implies that all the contributions stemming from 
integrations are subdominant in respect to the first two lines in (I4.62I ). this provided that \fI*_ {±q)\ < 1/2 for 
all configurations of variables in {z+j that belong to {Aq, ±q]. This condidtion is not satisfied, especially if |{z+)| 
becomes large. One should then invoke conjecture IB. II stating that, in fact, higher order oscillating terms in the 
representation (IB. 511 ) for //J,;^''^ are more dampen than it is apparent from the sum in (IB. 511 ). This is enough 
to show that, indeed, the contributions to the oscillating tems at -2pf, u{Ao) - u{q) and frequencies that are 
stemming from are subdominant in respect to the terms appearing in the first two lines of (14.621 ). There will 
of course arize terms oscillating at higher multiples of these frequencies exactly as it happens in (IB. 511 ). These 
higher oscillating terms give rise to other critical exponents. For instance, one can convince oneself that, from 
this structure, one recovers the whole expected tower of critical exponents for the terms corresponding purely to 
oscillations at integer multiples of u{q) - u (-q) as predicted in [2] on the basis of CFT-based technique. ■ 



Conclusion 

In this article, we have continued developing a new method allowing one to build two types of series of multiple 
integral representation for the correlation functions of integrable models starting from their form factor expansion. 
One of these series which we called the multidimensional Natte series yields a straightforward access to the 
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large-distance/long-time asymptotic behavior of the two-point functions. In this way, we were able to extract the 
long-time/large-distance asymptotic behavior of the reduced density matrix for the non-hnear Schrodinger model. 

In order to provide applications to physically pertinent cases, the method we have developed has to recourse 
to a few conjectures. The first one is relative to the convergence of the series of multiple integrals representing 
the correlators. This conjecture is supported by the free fermion case, where the convergence is rather quick, 
especially in the large-distance/long-time regime. The second conjecture concerns the possibility of using an 
effective series instead of the one appearing in the form factors expansion of two-point functions. Both series have 
been assumed to have the same thermodynamic limit N,L ^ +00. This conjecture is supported, on the physical 
ground, by the argument that sums over states whose energies scale as some power of the system-size ought to 
give a vanishing contribution to the sum over form factors once that the thermodynamic limit is taken. It would 
be very interesting and important from the conceptual point of view to prove these two conjectures in the case of 
models that are away from their free fermion points. 

However, we do insist that we have organized the analysis in such a way that all of the aforementioned 
convergence issues are separated from the asymptotic analysis part. Therefore, all the part of this work related 
purely to the asymptotic analysis is rigorous. Moreover, we do expect that the scheme of asymptotic analysis 
we have developed can be applied in full rigor to many cases which are free of convergence issues. We do also 
stress that, for the moment, the proofs of convergence of series of multiple integral representations for correlation 
functions of models away from their free-fermion point are, in general, an open problem. Apart from very specific 
representations related to the spin- 1/2 XXZ chain, the proof of convergence of a series representation for two-point 
functions could have been carried out only in the case of the Lee- Yang model by F. Smirnov. 

We have chosen to develop our method on the example of the one-particle reduced density matrix in the non- 
linear Schrodinger model. The case of the current-to-current correlation functions in this model will appear in 
ll65l . It seems however that the method is quite general and applicable to a vast class of integrable models where 
the form factors of local operators are known. In particular, it should be applicable not only to lattice models where 
the form factors admit determinant-like representations lfT4l l53l 1611 iTBll but also to integrable field theories where 
the form factors of local operators can be computed through the resolution of the so-called bootstrap program. 
For instance, the method seems applicable to the analysis of certain two-point functions (and their short-distance 
asymptotics) in the sine-Gordon model whose form factors have been obtained in f8ni821. In the latter case, we 
expect to deal with some multidimensional deformation of the 3'^'^ Painleve transcendent, a new type of special 
function whose description and asymptotic behavior is interesting in its own right. 
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A Thermodynamic limit of the Form Factors of conjugated fields 
A.l Thermodynamic limit of form factors 

It has been shown in f64l with the help of techniques introduced in H6l ITSl that the normalized modulus squared 
of the form factor of the conjugated field taken between the ground state {Aa}^ and any finite n particle/hole type 
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excited state {juij'^'^^ admits the below behavior in the thermodynamic limit N,L ^ +00, N/L D 
|(^(Kri)|cDt(0,0)U(M)f)>f 



1) 



Above, Fq corresponds to the thermodynamic limit of the shift function associated to the excited state of interest 
(at y6 = 0). The auxiliary parameters of Fq are undercurrent by the various functional appearing above. We recall 
that the parameters fik appearing in the rhs of (lA.ll) are defined as the unique solutions to ^ - a/L. 

The discreet part 

The first two functionals appearing in (lA.ll ) correspond to the leading in L behavior of the so-called singular part 
^N[F{c^],^{ej,^] of the form factor, namely 

Given any function v{A) holomorphic in some neighborhood of [ ; ^ ] , one has 

The parameter A^+i appearing above is defined as the unique solution to L^y (/Ia^+i) - N + \,k[v\{A) is given by 
(14.471 ) and G stands for the Barnes double Gamma function. Finally, we agree upon 



pr / sin [7rv(///,J] \^ FT ( " ^ " ^ + ^(j"p«)' P«, + 2 - /j^ - v C"/,J , ha + vijiij 

U\ TT / 'U I Pa-N-l, Pa + v(MpJ, N + 2-ha, ha 



ha -Ph 



0=1 ^ ' £1=1 

(A.4) 

There 

SM(.) = 2v(a.)ln(^^).2f^:^4^di and .^A,,) = 2n-p^. (A.5) 
\(p{aj,-q)) J A-oj p{A)-p{p) 
-1 

Above, we have used the standard hypergeometric-type representation for products of F-functions: 
ai, . . . ,a„ \ 



F 



n — • (A-6) 
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Description of Qf\ 

In order to give an explicit representation for Qf}y we need to introduce a few notations. First, let Fp correspond 
to the thermodynamic limit of the /3-deformed shift function associated to the choice of the rapidities {jJ-p^} for the 
particles and {/i/,^} for the holes. The auxiliary arguments of the shift function will be kept undercurrent. Also, let 
ni € N and Us be the open strip (12.151) around R. 

Then there exists yo > small enough and y6o £ C with % (/3o) > large enough and 3 (/3o) > small enough 
such that 

^2inyFp(o^)_y ^ q Uniformly inn - 0,...,m and [y,co,fiAyipXiMrX^ e D^^^xU sxU p.xU'lxU'l . (A.7) 

Let all parameters /i/,^, a = I, . . .,n belong to a compact Kq+^ D [-Q^q] for some 6 > and let 'i^q+f; be a small 
counterclockwise loop around this compact then admits the below representation 



[ f^h„ -q + eic e^i-c[yF,]{^„^^ac) | ^^^^^ 
I fip^^ - q + eic e2™C[y/7j](^,,„+eic) I 




^ - ^lh, - ic){^K - Hp, - ic) detc^... + YViyF/,]] detc^^,, [/ + rVjyF^]] ^^^^ 
}Ji iP^Pa - f^Pb - ic) (Ph„ - - ic) det^ [/ - Kjln} 

There C[F^] stands for the Cauchy transform on [ ; ^ ] and Co[Fyj] is given by a double integral 



CyFp] {A) = and Co = - — rd/ld// . (A.9) 

J 2m /u-A J n-u- ic)^ 



{A-fi- ic)' 



(^-^^pM^-^^K^ic) \ cr2,>.viM-cf2,-.vi(..,.) ^(^-^0 ... 

{o)- ^ihJ[oJ- Hp^+ic)^ e 1 

The representation (IA.8I ) is valid for « = 0,...,m and (7,^8, {//pj", ^ ^^0,70 >^ Ujs^ x x K^^^ and 

defines a holomorphic function of these parameters belonging to this set. 

It is also valid at 7 = 1 , provided that % (fio) > is taken large enough for condition (IA.7I ) to be fulfilled at 
7-1. 

Finally as follows from proposition lA. II given below, the product Do;L[F/3]'^A',n[F/?]6',*f | is holomorphic in 
5lyS > 0, and can thus be analytically continued from t/^g up to /? - 0. It is in this sense that the formula (I A. II) for 
the leading asymptotics in the size L of the form factors of is to be understood. 

Proposition A.l 

Let m e N, 5 > small enough define the width of the strip Us around R and {[jJLpXl'^ {^J■haYlj £ U'^ x K^^^, 
where e > and the compact Kq^^ is as defined by (14.271 ). 

Let V, h and t be holomorphic function in the strip Us around R and such that h{Us) c [z. : %{z) > 0] and 
3 {h (z)) is bounded on Ug. Set vp {A) = y (/I) + ifSh {A). 
Then, there exists 
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• ySo € C with % (J3q) > large enough and 9 (/So) > small enough 

• To > and small enough 

• a small loop ^q+e Ug around the compact Kq+g; 

such that uniformly in /3 € f/gg the function A h-> Q-^"^~i'('^+'l^f''>W _ i has no roots inside o/'^^+f. In addition, the 
function 

n 

{M\MX,,y,IS) ^ G(l -rT(-^))G(2 + yT(^))]~[(e-2-^^/^(^'0 - l) • det^,,, [/ + 7^ M ({yUpJ?, Kli )] 

a=\ 

(A. 12) 

is a holomorphic function in Ug X K'^^^^ X Do,7o x f/gg, this uniformly inO < n < m. 

It admits a (unique) analytic continuation to U'^ X K^+^ X i)o,i X |z e C : %.{z) > 0). In particular, it has a 
well defined /? ^ limit. The /? ^ limit of this analytic continuation is still holomorphic in ({^pj", [l^ha]" ) ^ 

In (IA.12I ) we have insisted explicitly on the dependence of the integral kernel "V on the auxiliary parameters 
{[ppj", {^hj" ), cf (lA.lOl ). The same proposition holds when the kernel is replaced by 'V as it has been defined 
in lTOT]) . " 



Alternative representation for Q^-- 



It so happens that the smooth part of the form factor's asymptotics admits a representation as a functional acting 
on a unique function H. More precisely, 

Lemma A.l Let m € N and the strip Us be fixed. Let A > be some constant defining the size of the compact Ka 
(I4.27I ). Then, there exists A, 5, m dependent parameters 

• jSo e C with % (J3q) > large enough and 3 (ySo) > small enough 

• To > small enough 

such that uniformly in ({ja)", {Za)") E K'^x K'^, \y\ < yo, /? € t/gg and < n < m 



OS) / M"] _ c^iP) 



{Za\\ 



{ya\\ 

[ZaXl 



with H\A 



{Zaf, 



" 1 

y— 



(A.13) 



The functional acts on a bounded loop 'tf (K^) Us around the compact K^. The functional [m] is 
a regular functional ( cf definition \4. 1 D of m in respect to the pair {Mc^^,'i^ (Ka)) where the compact Mc^^ has its 
boundaries given by 'rfgut cmd cis depicted in the rhs of Fig. [7] For all m € G (Mc^^) such that ||<:^(if^) < C^^, 
where Ci^^ is a constant of regularity of the functional one has 



det.^^ [/ + [yG^, m]] det<^^ [l + y~f \yGfi, xn^ 



X exp < - 



dydz 
{linf 



m{y)'uj (z) In (z - 3^ - ic) 



det^ [/ - Klln\ exp Vlin X C[yG/j] {q + iec) 



(A. 14) 
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In the above formula, one should understand Gp as a one parameter family of functionals of vj given by 

Gp{A) = Gp{TnMA) = {il3-\l2)Z{A)-cp{A,q)- ^vj{z)(p{A,z) . (A.15) 

In the second line of (IA.25I) there appear Fredholm determinants of integral operators acting on a contour 
The contour corresponds to a loop around ^ {Ka) such that ^ c J/^. The kernels read 

nv,ru^{<.,c.') = Zi^:izl_exp( (f (z) In ( ^^^] \ • e^[^'--JM-C[2,-.v](...) -^0 ^^^^^ 
2n a> - q + ic [J 2in - z. + ic I j Q-j.mv(aj) _ i 

and 

nv,^M<.,u') = J-^^^expj (fi i^^(,)ln(^^i^)|.e™-)-™--^)^|g^ (A.17) 
2k CO - q - ic [J 2in xo) - z - icl } Qnmco) _ i 



The A, m and 6-dependent parameters /3o and yo and the compacts 'if (Ka), Mc^^ are such that the constant of 
regularity C's^ satisfies to the estimates given in (14.301) and is such that one has 



V \\Tu\\^iKA) < C'Sa ||roG/J < 1/2 and < Os^ uniformly in ({3;,}'^, {zaT^ e x ■ 

Proof — 

We first clieck tliat is a regular functional. 

• Gyj [ct] is a regular functional as it is linear in m and {Ka) is compact. 

• the estimates le'^ - e-^'l < e'"^'^'- ' \x - y\, majorations of integrals in terms of sup norm and derivation under the 
integral sign theorems ensure that all of the exponential pre-factors in (IA.25I) are also regular functionals of 

TU. 

The associated constants of regularity can be taken as large as desired. It thus remains to focus on the Fredholm 
determinants. For this let us first assume that we are able to pick the contours '^outjin delimiting the boundary of 
the compact M<^^ in such a way that there exists 

ySo e C and yo > such that Q^i^yGA^(*,y)\{X) _ 1 ^ q V (i, j,yS, y) e t/^ x x 11^^ x 2)o,7o (^-18) 

this for any holomorphic function w{X,y) on Mt-g^ x W^, c C^' , that satisfies ||tn'||<^(if^)x»'^, < C^^. 

If this condition is satisfied, then the integral kernels [yG^g, yG^, a>') and yi^ [yG/3, yG^, vj\^(X), oj') 
are holomorphic in o), cd' belonging to a small neighborhood of and y € Wy. The contour 't^a being compact, 
the two integral operators yV [yG/j, yGyj, tirj and yV [yGyg, yG/j, cij are trace class operators that have an analytic 
dependence on j € Wy. 

Recall that if A, B are trace class operators (|Hli stands for the trace class norm) then 

|det [/ + A] - det [/ + B]\ < \\A - B\\i e"^"'+"^"'+' . (A.19) 

Also f76l, if A(j), y e Wy c C^>, is an analytic trace class operator then det [/ + A(y)] is holomorphic on Wy 
These two properties show that, indeed, in (IA.25I ). the two Fredholm determinants of integral operators acting on 
the contour are regular functionals of m. 
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Hence, it remains to prove the existence of yo and such that condition (IA.18I) holds. Given m{A,y) € 
ff{M<^^ X Wy), the function ca ^ Q2i7TyG/,ivT](^) _ ^ ^as no zeroes provided that 

\yoGfi [m{*,y)] {A)\ < 1/2 and 3 (G/j [m{*,y)] (A)) > uniformly in {A,y,/3) e U^xWyXU/}, . (A.20) 

One has that, for /3 e U/jg 

5(G^[m](A)) >%(fio) inf [%{Z(A))]-id (fio) + l/2)p {Z)\\ -\\cp\y_-\\rjj\y^^^^^ sup . 

Hence, 3 {pn \w\ (i)) > as soon as ll^cr||>^(jf^)xvv'. — Cc^^ with 

C^, = sup ^ ^ 10 z)\ I • OSo) mf^ [51 (Z (A))] - (0 OSq) + 1/2) ||3 (Z)||^^ - \m^2^ . (A.21) 

Here 21 (/3o) > is taken large enough for C^^ as defined above to be positive. Then, if ||ti''||'^(ifA)xH'j, ^ with 
Cff^ as given above, one has 

sup \Gp[m(*,y)] (aj)\ < (1021 OSq) + 3 CSq) + 1/2)||Z||[;^ + ||<^||j;2 + Wmiy^KM sup (£ ^ \(P(A,z)\ 

coeUs " AeUs J 

yeWy nKA) 

< (imOSo) + 23 0So)+l)l|Z||f;, . (A.22) 

Hence, if we take y'^ = 2(1121 (/3o) + 23 (ySo) + 1) \\Z\\u,, the condition \jGp [m]\ < 1/2 will be satisfied for all 
lyl < To and p e t/^^. It remains to tune 21 (fio) so that conditions 

7rd(5M%,^(^A)) 2m 

C% • 1 r^^ — — > A and < C<^, . (A.23) 

^ \dM^^\ + 2ndidM^^,'^iKA)) d{KA,'^{KA)) 

are satisfied. 

One can always choose the contours %jut/m defining 5M^^ in such a way that d(5M^^, ^ (^a)) > c this 
uniformly in A>0. These contours can also be chosen such that there exists an A-independent constant ci with 
I^M^^I < ciA. It is also clear that the contour '^{Ka) surrounding the compact Ka can be chosen such that 
I'l^ {Ka)\ < C2A for some A-independent constant C2 and also d{KA,'rf (Ka)) > c' . It is then enough to take 
2^0So) > C/JqA^ with c/Jq being properly tuned in terms of c,c' ,c\,C2 so that conditions (IA.23I ) hold for any A 
sufficiently large. 

Note that the second condition in (IA.23I ) guarantees that the function H as given in (IA.13I ) satisfies < 
C<Sa uniformly in the parameters ({yUpJ", ljU/,J'^) ^K\x.K\. 

Having proved that is a regular functional with a regularity constant C<g^ > sufficiently large, we can 
evaluate it on H. Then, it is readily seen that Gyj [m] (A) coincides with the shift function F/^ once upon taking 
m = H as given in (IA.13I ). All offier integrals involving m - H we. computed by the residues at //p^ and All 
calculations done, one recovers ffie representation (IA.8I ) for the function Q^^y. We stress that the parameters yo and 
y6o ensuring the regularity of the functional are also such that ^J^], is well defined due to conditions (IA.20I ). ■ 
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Regular functional for Q^-y 

A very similar representation to the one given in the previous lemma exists for the functional 0N-y 

Lemma A.2 Let m € N and the strip Us be fixed. Let A > be some constant defining the size of the compact 
(I4.27I ). Then, there exists A, m and 6-dependent constants 

• ySo £ C with % (J3q) > large enough and 3 (y6o) > small enough, 

• To > small enough, 

such that for L large enough and uniformly in [iHpJ", {/^/i^}") € K^X K^, \y\ < yo <^nd < n < m 



The functional "^^^ acts on a bounded loop {Ka) c Us around the compact Ka- The functional [cr] is 
a regular functional ( cf definition \4. 1 D of m in respect to the pair (Mc^^, ^ (A'^)) where the compact Mcg^ has its 
boundaries given by '^out cmd ci^ depicted in the rhs of Fig. [7] For all m ^ G (Mc^^) such that ||w||<^(/f^) < C^^, 
where C-^^ is a constant of regularity of the functional ^^1^' 

M = Wiv[rG^] ( ) n I - (z) {- In (^iV;.[rG/j]) (z) + In (z - Pn^i + iec) 

( fdvdz ^1 det% [/ + 7'^[7G/},rG/},TZj]]det.^^ I + yf N[yGp,yGp,m\ 

xexp<- (f) ' m (y) Tu {z)ln{z - y - ic) 



(2^^)' J detA,+i[5(A')[^]]det^[5W[^^Gj] n V-.^rG^[Ti7]]0u;v+i) 

(A.25) 



In the above formula, one should understand Gp as the one-parameters family of regular functionals of m as 
defined by (IA.15I) . We did not make the functional dependence ofGp on m explicit in (IA.25I ). The functionals 
and Vn;(; have been defined in (I4.6I ). We have added the [yG/s] symbol so as to make it clear that the parameters 
[Aa]^ entering in their definition are functionals ofyGp through the relation Aa = (a/L). 

In the second line of (IA.25I ) there appear Fredholm determinants of integral operators acting on a contour 
^A- The contour '^a corresponds to a loop around ^ {Ka) such that ^ c Ug. The kernels read 

In CO - fiN+i + ic [ J 2m \a)-z + ic/} VNo[v]{a)) e -^"'y('^) - \ 

nKA) 

and 

In CO - - ic [ J 2in \co-z-ic/} VNo[v]{to) Q^^imoj) _ { 

The constant of regularity C^^ satisfies to the estimates already given in (14.301 ) and is such that 

V \\m\\^(KA) < C^A \\yoGp Mlljy, < 1/2 and < C^, uniformly in {{/UpJ",, {phj") ^KlxK^. 
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Proof — 

The proof is very similar to the one of lemma lATT] Hence, we only specify that for L-large enough, and as soon 
as condition |yG^ [cj] {A)\ < 1/2 for all A e Us is satisfied, the parameters Aj are seen to be regular functionals of 
vj thanks to their integral representation 



All other details are left to the reader. 



A.2 Specific values of the functionals ^^^j and 

08) 



In this subsection, we estimate the value of the functional '^y.^ [tc] for a specific type of function m. This result 
will play a role later on. 

Lemma A.3 Let the function v{A) = v{A | Izk)" Ayk]"^^ ) be the unique solution to the linear integral equation 
driven by the resolvent R of the Lieb kernel (ie (/ - K/2n) (I + R/ln) - I): 

y(A) + y I ^R{A,^x)v{p)^{ip-\|l)Z{A) + J]cf>(A,Zk) - J] ^''t) ■ (A.29) 

Let A > be large enough and such that {^{Zk}" , {j/t)"^') £ X ^a^'- Po ^ C andjQ be the two numbers 
associated to the constant A as stated in lemma \AA\ Then defining 

n+l 

^ = y, T— - > : T— - I -^-^dr , (A.30) 



j^^A-Ja A-q ^^A-Za J {A-Tf ^' 



the below identity holds 

n n+l P 

n n (yt-Za- ic) iza -yb- ic) det„ \6k( + jVkc [yv]\ det„ 5ke + yVke [yv] 



[^] = -ic 



n (ya-yb- ic) n (Za -zb- ic) 

a,b=\ a,b=\ 



(A.31) 



The non-trivial entries of the two determinants are given by (14.81 ). The auxiliary variables {{Zk}'l, {yk\'l^^) on which 
these entries depend are undercurrent by the set of auxiliary variables on which depends v. 

Proof — 

The function v is bounded on the strip Ug- As a consequence, the associated function m (IA.30I ) is also bounded 
by an A independent constant. The estimates (14.301 ) for the constant of regularity Cc^^ for the functional ensure 

that there exists A large enough such that ||?jt||<^a:^) < C^^ uniformly in {{zk]"^, {yk]'l^^) e K^xK']^^. Thus, '^^l [m] 
is then well definied. 

A direct calculation leads to 

n+l 



exp 



zj{z)\n( "^"^ )| = n "^"^^ n ^^i^i^eC[2-H(-±-)-c[2,>HM . (A.32) 

2in \(JL> - z±icl to - q ^\a)-ya±ic'^l oo - Za 
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By using the linear integral equation satisfied by v and the representation (IA.15I) we get that 

Gp\m]{X) = v(.A). (A.33) 
As a consequence, the kernel Y and "V simplify 

T yGfl [ct] , TET (oj, a> ) = - ; — - — (A.34) 



and 



T \yGR[m\,m\{(jj,(ji)) - \ —■■ — ■r-^ r (A.35) 



The associated Fredholm determinants can now be reduced to finite-size determinants by computing the poles 
at 0) = Za with a - \,...,n (by definition of 70 and ySo, since I7I < yo and yS € U/j^, there are no poles of e^"^''^'^) - 1 
inside of '^a)- 

This leads to 

det% [7 + rnyGf} [m] , ru]] = det„ [Ste + yVu [yv] [{Za}\, {yaW')] (A.36) 
det% [7 + yr[yGp [tu] , m]] = det„ [fe + yVki [yv] {{Za}", {ya}T')] ■ (A.37) 

The claim then follows once upon applying the identity 

n n+1 

n U(yb-Za- ic) (Za -Jb-ic) ( rA A 

- ic — ; ^ exp - (p -m (y) m (z) In (z - j - ic) 

n ija -Jb- IC) n (Za -Zh- IC) ^eiKA) 
a,b=\ a,b=\ 

X eCo[7G/,[^]]g-2'^ 2^c[yG;,[t^]](?+<6c) | _ (^^'uj (z) {2ctC [tG^M] (z + /ec) + In (z - ^ + /fc)} 

(A.38) 



Lemma A. 4 Let y be small enough and L large enough such that v^^^ (jj.) is the unique solution to the non-linear 
integral equation 



V 



N+\ n+1 N 

a), 



\A) = {il5-\l2)Z{A) - c^{A,q) + Y,(^{A,^ia) - Yj<P^^^ya) - (A.39) 



a=l a=\ 0=1 



The parameters Aa appearing above are functional of v^^-' through the relation ^yy(L)(Aa) = a/L, fia are such that 
^ (Ma) = cijL and the parameters ya e t/^/i cire arbitrary. Finally, L is assumed large enough so that all parameters 
[^a]^, {3'a)i^') of the IN + 2 + n-uple belong to K2Al- Then, given /3o andyo as in lemma \A2\ one has 
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the identity 



HI* 



1 



-ic- 



n n+l . 

nn(yb- K 

a=l b=\ 



ic) {^i, -yh- ic) 



n+l 

n 

a,h=l 



i+l n . . 

n {ya -yh- ic) n [K - K - 



a,h=\ 



X det„ [6u + yVkAjy^"^^] (UiXv ^ya)T')] det. [he + yV^drv^^^] (UiJv Mr') ] (A-40) 



Proof — 

It has been shown in proposition ID.3I that for lyl small enough and L large enough the solution v'-^^ to the 
non-linear integral equation occuning in the rhs of ( IA.39I ) is unique and exists. Moreover this solution is bounded 
on Us by an L-independent constant. 

As discussed in the proof of lemma IATTI the contour '^(^2Al) can always be taken such that, uniformly in L, 
d{^{K2Ai) > Kxal) > c' > for some constant c'. Hence, the principal argument A of His uniformly away from the 
compact Kial where the auxiliary arguments of H are located. As a consequence, it follows from the expression 
for H and the estimates for the spacing between the parameters /Iq and Aa 



fia-Aa = Inyv^^^ ifia) I (Lp'ifi,)) + O (L"2) , uniformly in a - 1, . . . , (A.41) 

that ||^||<i^(/f2A ) bounded by an L-independent constant, this uniformly in L large enough. In particular, for L 
large enough, due to the estimates (14.301 ) for the constant Cg;,^^ of regularity for we get that ||^||^(^^) < 

C^2Ai • One can thus acts with the functional ^^^^ on H. A straightforward residue calculation shows that 



HI* 



N+l 



(A) . (A.42) 



This means that all the Aa appearing in the expression (IA.25I ) for the functional ^^-Ial coincide with the 
parameters Aa defined above. The claim of the lemma then follows from straightforward residue computations 
and multiple cancelations. The Fredholm determinants reduce to finite rank determinants that can be computed 
by the residues at w = Ai^, a = I, . . . ,n. ■ 



A.3 Leading asymptotic behavior of one particle/one hole form factors 

We now build on the formulae for the leading asymptotic behavior of form factors so as to provide, properly 
normalized in the size of the model, expressions for the large-L limit of the form factors of the fields between 
the A/^-particle ground state and N + 1 -particle excited states corresponding to one hole at one of the ends of the 
Fermi zone and one particle either at the other end of the Fermi zone or at the saddle-point Aq of the function 
u (A) given in (12.141 ). Such thermodynamic limits of properly normalized form factors appear as amplitudes in 
the large-distance/long-time asymptotic expansion of the reduced density matrix. The exphcit expressions that we 
write down will allow for such an identification. We do stress that all shift functions appearing below are taken at 
y6 = 0. The fact that (IA.46I )- (IA.48I ) are well-defined in this limit follows from proposition lA.ll 

In the following, let {A} = [Aa]^ stand for the Bethe roots corresponding to the ground state in the A/^-particle 
sector. Let {yu®} = {yu®}^''"^ stand for the Bethe roots corresponding to the ground state in the (N + l)-particle sector 
Taking into account that F'^ stands for the thermodynamic limit of the corresponding shift function cf (13. 51 ). we 
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define 



N,L-*+oo\27T ) 



(«A({/^^|)|ot(0,0)|^(|i|)> 



(A.43) 



Similarly, given the set {yu^^} = ilJ-q'^}^^^ corresponding to a particle-hole excitation such that pi = and hi = 



N + l,v/e denote by Fq the thermodynamic limit of the corresponding shift function cf (13.51 ). and define 



(«A(|y"-''|)|cDt(0,0)U(W)> 



(A.44) 



Finally, given the set {ij.q°} = {jU^")^^^ corresponding to a particle-hole excitation such that hi = N + I and 
f^Pa - ^0 we denote by Fq° the thermodynamic limit of the corresponding shift function cf (13.51) . and define 



IT' 



lim — 

N,L~*+oo\27:/ 





4°])\^H0,0)\^{{A])) 




<A(|il) 









(A.45) 



By using (lA.ll) and expressions (IA.3I) - (IA.5I ) we are lead to 

|r,-t = [f;'. p] s [f;-. p] gf, ( -« j exp {,| ([f;«(-,) - 1]^ - [f.-'w j^)) . 

and finally 

= K> /'l ^ K P] 4;! ( ) -p {/^ in-,) f - [Fti,) + 1]^)} . 
The functionals S, J{± and appearing above have been defined in (14.461 ). (14.471) and (I4.48I ). 



(A.46) 



(A.47) 



(A.48) 



B The generalized free-fermion summation formulae 

In this appendix, we establish summation identities allowing one to recast the form factor expansion of an analogue 
of the field/conjugated-field two-point function that would appear in a generalized free fermion model in terms 
of a finite-size determinant. The representation we obtain constitutes the very cornerstone for deriving various 
representations for the correlation functions in the interacting case. In particular, it allows one for an analysis 
of their asymptotic behavior in the large-distance/long-time regime. We first establish re-summation formulae 
allowing one to estimate discreet analogs of singular integrals. This will open the way for obtaining Fredholm 
determinant like representations out of the form factor based expansions. 
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B.l Computation of singular sums 

Let ^ stand for the thermodynamic Umit of the counting function (12.61 ) and E- be a non- vanishing and holomorphic 
function in some open neighborhood Us (cf (12.151 )) of R such that 51 (in EZ^^ has, at most, polynomial growth, ie 

1^1 [in EZ^ (A) - /Ci/]| < C2 1^1 (//"^)|+C3 , for some Ci , C2 , C3 e R+ and A: > 1 uniformly inAeUs ■ 

(B.l) 

We remind that the neighborhood Ug is always taken such that ^ is a biholomorphism on Ug 
In the following, we study the below singular sums over the set {//q}: 

<S^^^[£:^1 (A) ^ V , TT7 with fia being the unique solution to ^Qia) = a/L . (B.2) 

^;^^2;rZ^'0u„)(//„-i) 

The summation runs through the set Sl = {a eZ : -wl < a < wl) where wl is some L-dependent sequence in 
N such that L - o (wl) and [wl • ^ I =0 (L). 

Proposition B.l Lef A^^, be a compact neighborhood of[-q;q] lying in Us, then under the above assumptions 
and provided that L is large enough, one has, uniformly in A e Nq 

Sf[E-.'YA) = ^f^EZ\p) + lf\E-_^YA) (B.3) 

'^bk;L 

.nr 71 5 rduE-^iu) SaIez^A)] E-\A)L^'(A) a)t 71 

S\^>\E-^{A) = — Jt^^ _ I L — J + ;r " + /r (i) . (B.5) 

2 L - ^ 5i j 27r //-^ e2-^('i)-l 2sin2[;rZ^(i)] 2 L - ^ 

r/ie integration goes along the curve %,k\L depicted on Fig. |?] Also, given r e N, 

,(L)r^-2i.,._ Cdz Ez\z) 1 r dz g:^ (z) 1 f dz £:^(z) 

^ -J^^ j 27r(z-i)'-l -e-2'>^^« j 27r(z-i)'"e2'^^(J)- 1 j 2n{z-Af ^ ' 



^T;L '^ii 



L 



r/je contours '^yi-i cifs depicted in Fig. \3\whereas '^hd;L is depicted on Fig. |?] 

The junctionals tP\^z}\ {A) are such that lf^\^z}\ {A) = O {(L/wl)^^'''^), uniformly in A e Ng. 

Proof — 

Let Nq be a compact neighborhood of [ ; ^ ] in Ug. Then, for L large enough it is contained inside of the 
contour 't^^-l U 'ifi-i as depicted in Fig.[3l and thus 
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Figure 3: Contour 'ia-^-^i U ^j^;/, lying in Us. 




Figure 4: Contours '^hk;L (solid lines) and 'lahd^i (dashed lines) in the case of k odd {Ihs) and k even {rhs) both in 
the case C\ < 0. The dashed lines 'lohd-L are pre-images of the segments [ eyWi + eyl2 ; €yWi + e^ll + ie'a], with 
V € |1, r) and ei = -1 and = 1. The sign of e' depends on the left or right boundary, the parity of k and the sign 
of Ci. 



iEz\A) r EzHz) 1 

dz 



Ez^iz) 1 r £"1^(2) f e2'^^(^) ) /£':^(^) 

- U dz - 



r Ez\z) 1 r £:2(z) f e^'^^^^) ) 

J 2n{z- A) e^inmz) _ i ^ + J 2n{z-A) [e^inmz) _ i j ' 



-2 



d;u£:^(/^) . EZ^iA) 



Jin IX -A 'e^inmi) 

^bk;L 



^ + lf\EZ^]{A) . (B.7) 



In order to obtain the last line, we have deformed the contour '^■^■l into the contour 'i^hk;L U '^hd-L as depicted in 
Fig.m The intermediate points ±A entering in the definition of '^m;L are chosen large (in order to include Nq) but 
fixed, in the sense that L independent. 

The representation for •S2^^[£'Z^ j (A) follows by differentiation. The computations for iSg'^^^f'Z^ j (A) are carried 
out similarly with the sole difference that there is no pole at z = /^. 

In now remains to prove the statement relative to the asymptotic behavior in L of the functionals /J-^^I^SZ^j. 
The main difficulty is that the function E- (A) might have an exponential increase when A belongs to the upper or 
lower half-plane. We establish the claimed estimates for the '^-|-;i,-part of the contour. This can be done similarly 
for 't^i-l, and we leave these details to the reader. 
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We first perform the change of variable^ z - ^'^(s) and set ul - wl/L + 1/(2L). The contour of integration 
is then mapped to the contour depicted on the rhs of Fig. [3l We stress that the parameter a > is chosen in such 
a way that "W^-l U '^ul lies in Ug. The aforementioned change of variables leads to 



r dz EZ^ (Z) 1 _ f ^ ^ ^-1 r • ^ 1 _ .ni^a - 



a 



4 o Z7-^ 1 

^ I _ ^-,2nLa ^—2msL\ 

[r 1 + ia) - A^' ^ ' 

'^^orH/. + .«L)--— ^— — 7{l+e2-^)-^ . (B.8) 



We first establish a bound for the integral over the line [ -ul ; ml ]■ It follows from the integral equation (12.71 ) 
satisfied by p that, 

p{X) = X±nD- IcDIA + O (^"^) when 51 (i) ^ ±oo . (B.9) 
Hence, uniformly in < t < a and for 5 € R large, 

^'\s + h) = >f/s + 2ott + 0{s~^) where i/^,. ^ 2;:^ - ;rD(l ± 1) + — € R . (B.IO) 



The condition (IB. II ) implies that there exists constants C > 0, C > such that 

\% [ln£:2(/l)]| < C |5l (//)| + C , uniformly in i e [/^ . (B.ll) 

As a consequence, uniformly in < t < a, 

\^ [inS:^ o (s + /t)]| < C |3 [(iA.v)*^ + 2inkT{il/sf'^ + 0(/-2)]| + C 

< CTA:(27r/|5|*^-^|3[l + 0(5"^)]| + C' . (B.12) 

There exists an sq such that |0(5~^)| < 1 for 1^1 > ^o, this uniformly in < t < a. Moreover, for such an sq, 
we define 

C" = C + max|5l[ln£':2o^-i(s + /r)]| , (B.13) 
with the maximum being taken over 1^1 < sq and < t < a. Hence, for any s eM. and < t < a 

|2l [inS:^ o {s + /t)]| < 2kCa{2nf Uf-^ + C" . (B.14) 
Therefore, we obtain the estimate 



As E-^ 



J 



iTiLa -linsL 



2n f [^-1 + 'a) - /I] 



e^" ] 2wL + l e2*^c°^(2;r)*„; 



k„k-\ 



' we remind that <f is a biholomorpliism on and that p' > of 1 
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where we have used that {wl ■ L ^)'^ ' = o (L). It remains to estimate the integral over the Unes [ ; a ] : 

v-i r i'^sEJ I 1 

Zj^ ^ — 0^ + f"L)7;rT-^ 

J In ^' {^~^{is + e 



1 + e 



euL) - X) 

I \ 

sup < I 

.ve[0;a] ^' o ^-1 (/^ + euA [^"1 (/5 + eM/,) 

ee{±l) 



^2Ck{2nf{uLf-'r 



+ e 



2;rLT 



-dr . (B.16) 



By making the change of variables y = Lt and then applying Lebesgue's dominated convergence theorem, one 
can convince oneself that the integral in the second line of (IB. 161 ) is a 0(L"^). 

The last class of integrals to consider stems from integrations along 'iobd-L- In order to carry the estimates, we 
need to use the finer condition (IB. II ). Here, we only treat the case of k even and C\ < 0. All other cases are treated 
very similarly. An analogous reasoning to (IB. 121) leads, uniformly in < t < a to 



2l[ln£:2o^-'(5±/T)] -T(+A;(2;r/Ci/-^ +0(/-2)) for +s>0. 



(B.17) 



There exists s'^ such that for 1^1 > s'^ one has |o {^s^ < k{2ni^ |Ci/ ^| /2. As a consequence, for \s\ > s'q and 
+s>0 



%[lnEZ^ o^~\s±iT)]<-k—^T\Cis''-'^\ uniformly 0<T<a. 
Therefore, 

f dr (^:'/r)°rH6ML-/eT) 



(B.18) 



/ 



In - nY 



< 2 sup 

re[0;ff] 
ee|±l| 



\^-\euL - ier) - a]" 



j £e-*^(2./"i"l^'l5 =0(4"^-'-) . (B.19) 



B.2 The generating function: form factor-like representation 

From now on, we assume that the function takes the form Ez'^ (A) = e'™('')+5('*) where u (A) is given by (12.141) 
and g is a bounded holomorphic function on the strip Us around R (12.151 ). We also assume that v € t& (Ug). 
We remind that the parameters {fJ.a}aez (resp. {Aa}aez) defined as the unique solutions to L^Qia) = a, (resp. 
L^v (Aa) = a), where ^ is given by (12.61) and (A) - ^{A) + v{A) /L. We define the functional [v, E^^ as 

N 

N+\ n E_ (Aa) . . 

Z Z 5„ ^^j; . <B.20) 

„=o pi<-<p„ ht<-<h„ r\ E^(uf) 



The functional Di^^„ has been introduced in (14.41 ). The sums in (IB. 201 ) run through ordered n - uples of integers 
P\ < • • • < Pn belonging to S^"' - N + l'^ and through ordered n - uples of integers h\ < ■ ■ ■ < K 
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belonging to = |[ 1 ; A'^ + 1 ]]. Finally, = |j £ Z : -wl < j < wl) and the sequence wl ~ . In particular, 
when L — > +cx3, wi grows much faster then N. The integers {pa] and {ha\ define the sequence {\ < • • • < ^at+i as 
explained in (I2.3I ). 

The functional Xi^[v, E^] admits two different representations. On the one hand, as written in (IB.20I ). X/v[v, E^] 
is closely related to a form factor expansion of certain two-point functions in generalized free-fermion models. 
On the other hand, after some standard manipulations |60|, one can also recast Xyv[v, f"^] in terms of a finite-size 
determinant which goes to a Fredholm minor in the N,L ^ +00 hmit. 

We derive this finite-size determinant representation for Xff[v,EV[ below. 



Proposition B.2 Under the aforestated assumptions concerning the functions E- and v, the functional Xj^\y,E^\ 
admits a finite-size determinant representation 



he + 



^'Mt) ^'Mt) 



(B.21) 



where 



P^^\A,h) = 4 
Also, we have set 



sin [ttv (i)] sin [nv (//)] 

2in {A - fi) 
sin [ttv (/})] sin [ttv (ji)] 
2^r 



£■_ (//) £■_ (A) ■ [0^^^[v, EZ^] (A) - 0^^^[y, EZ^] (p)] , (B.22) 
E. (A) iju) ■ 0^^^[v, EZ^]{A) ■ 0^^^[v, EZ^] (p) . (B.23) 



0^'\y,EZ'](A) = ip, 



'^EZStl + ^^^^ 

In p. - A 



-2mv{A) _ I 



+ 



(B.24) 



The contour of integration has been depicted on Fig. \5\and S^^^ (resp. 1^^^) is given by (Ib31) (resp. (IB^ ). 

'^^^■•^ - m-AL) = -WL- 1/2 



-3q 



mBL) = WL + 1/2 



3q 



Figure 5: Contour ^m;L appearing in the definition of 0^^\v,EJ-\{A), contour (solid line) and contour 
(dashed line). The contour ^m;L is such that, for > 4^ , it stays uniformly away from the real axis. 



Proof — We first recast the sum over the integers and {ha) corresponding to particle-hole like excitations into 
the equivalent sum over all possible choices of integers {a- < • • • < ^N+i with €a ^ = U S^"'. cf (12.31) . 
As all the sums are finite, there is no problem in permuting the orders of summation. Therefore, 

Y\El{Aa) 

^^[^'^-]- Z . (B.25) 
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The determinant entering in the definition of Djv can be represented as 



N 

n 

a=l 



■ detA 



1 



= det 



(1 - Sb,N+l) 
- -^b 



N+l 



da 



det AT 

|a=0 



1 



a 



(B.26) 



There we have used that for any polynomial Q of degree 1, one has 2(1) = 2(0) + 2'(0). 

It follows from the above representation that the summand in (IB.25I ) is a symmetric function of the N + \ 
summation variables /if^ that is moreover vanishing whenever = , k t a . Therefore, we can replace the 
summation over the fundamental simplex < ■ ■ ■ < ff^+i in the (N + 1)"^ power Cartesian product S^^^ by a 
summation over the whole space S^'*'^, provided that we divide the result by {N + 1)!. Once that the summation 
domain is symmetric, we can invoke the antisymmetry of the determinant so as to replace one of the Cauchy 
determinants by (A'^ + 1)! times the product of its diagonal entries. This last operation produces a separation of 
variables |[60l . Eventually, the result can be recast in the form of a single N x N determinant: 



4sin2 [7Tv{Aa)] 



with 



a=l ^vUa) 



■ y ^~ f 1 + — ) detA, \Mik + aPik (Un)] 



tS|^^^ being given by (IB.2I ). (IB.3I) - (IB.5I) and Pjk (jin) being a /i„-dependent rank 1 matrix: 

El (Ak) Sf\E-J] (Ak) 



InL 



(B.27) 



(B.28) 



(B.29) 



Using the fact that Pj^ (ju,,) is a rank one matrix that contains all the dependence of the determinant on the 
summation variable it is readily seen that 



(B.30) 



where 



EljAk) 
InL 



Sf^[E-_']iAt)-S\'^>[E-J]iAj) 



Applying (IB.4I) . (IB.5I ) and then using that L^{Ak) = L^y{Ak) - v{Ak) = k - v{Ak), we obtain that 



E^ {At) 



4sin^[7ry (/l^)] 



+ E.{^t)E-{h) 



linL {Ak - At) 



(B.31) 



(B.32) 



where is given by (IB .241) . Note that we have slightly deformed the form of the contours '^tk-L in 

respect to Fig. ID Very similarly, we find 



Jk- 



E-{Aj) E-[Aj)E-{Ak) 



0^'\v,E~J]iA,)0''^[v,E-J](Aj) = 



P^'\Aj,A,) 



'E-{Ak) 2nL L' J'"' L' J V AL^mnv[A^^mnv{Ak) 

where P^^^ {^Aj, Af^ is given by (IB. 231 ). It then remains to factor out the pre-factors from the determinant. 



(B.33) 
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B.3 Thermodynamic limit of [v, Ej^ 

Proposition B.3 The thermodynamic limit ofX^ [v, S^j is well defined and can be expressed in terms of a Fred- 
holm determinant minor. Namely, xAv,El\ X,A^Av,El\with 

q 

X^(c„)[y,£!] - <S^m[£:^] + 2 J ^sin2[7rvU)]F+U)£_a)0<^(„,[y,£:2](i) • det [/ + V] [y, £1^] . (B.34) 



1 ) 

2 , 



Here I + Y is an integral operator on \ \ acting on L ([ ; ^ ]) with a kernel 
sin [ttv (/I)] sin \nv (//)] 



V(i,/^) = 4- 



U_ {X) E. in) ■ {o^(oc) [y, (A) - [y, (^)| (b .35) 



2ot {A - ju) 

and the contour 'io''^^ dependent fiinctionals 0,^m) [v, j (/}) anJ "Sc^cu) [^^I^] ^ given by 



^ g-2(>y(i) _ 



Y and 5,^,H,[£:2]- . (B.36) 



(B.37) 



(/l) is the unique solution to the integral equatioi^ 

q 

sin [nv {A)] {A) + Jv {A, n) sin [nv (/i)] F+ {p) Ap = sin \nv {A)] E. {A) O^^^^y [v, EZ^] (A) . 

-'1 

Also, ^^"^^ - ^^""^ n [z € C : \'^iz)\<w]and ^^'"^ /jave Z^een depicted on Fig. |2] 

This representation can be seen as a generalization of the results obtained in i60ll . Also, the contour can 
be thought of as the L ^ +00 limit of the contour '^^hk■,L■ 

Proof — 

It is a direct consequence of the estimates obtained in appendix IB. II for iI.^\ez^] together with the fact that 
det^r ^dke + o (^~^)] ^ 1 in the case of remainders o(L~^) that are uniform in the entries, that 



Zw[y,£!l (s,^JeZ^] + -P\ ■ det[/ + y + aP] 



with I + V + aP acting on [ ; ^ ] and 

P (A, p) = - sin [Try (^)] sin [nv {p)] E. {A) 0") [y, EZ^] {A) [y, EZ^] Qi) 



(B.38) 



(B.39) 



Note that there is no problem with the integration over an infinite contour in 0,^(oo) ^v, E-^^ (y") and S^(oo) [e_^^ 

in as much as is built precisely in such a way to ensure the exponential decay of the integrand at infinity. 
Using that P is a one dimensional projector, we get that 

q 

det[I + V + aP]^det[I + V] \+a \iI + Vr\A,p)Pip,A)dAdp . (B.40) 



1+aJ (I + Vy'{A,p)P(p,A)dAdp 

^ -1 

It then remains to take the a-derivative and use the definition of F+ (A). 



' By no means F+ ought to be confused with the shift function 
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B.4 An algebraic representation for the Fredholm minor 

Proposition B.4 For L large enough, the finite N Fredholm minor [v, j defined in (IB.21I ) can be represented, 
through purely algebraic manipulations, as the below finite sum: 



n—\j ([,...,(;, 'yr(L) K—i a—i 

4sin2 [nvi\)] 



n 



yn+l - Zk 



I On+1 -^ik)(yk-zk) 



■ det„ 



1 



n 



(B.41) 



Above, appear two contours, which stands for a small counterclockwise loop around [-g',g] as depicted on 
Fig. \5\and "^^^^ = 'i^hk;L U ^T;i ^'^l;L^ 'Tobd;L U '^q- Note that '^bk;L it has been depicted on Figs. EH As shown 
on Fig. |5] '^^q stands for a small counterclockwise loop encircling 'tfq. Finally, the function f ^\y, v) is supported 
on 'tf^^'' and reads 

where 1a stands for the indicator function of A. 
Proof — 

The functional 0'-^^[v, EZ^] (z) as defined in (IB .241 ) is holomorphic in some sufficiently small open neighbor- 
hood of [ ; ^ ] . Hence, there exists a small counterclockwise loop around [-q',q] (cf Fig. |5ll such that the 
kernel V'-^\A,iu) admits the integral representation 



V'-^\A, p)^4 sin [nv (A)] sin [nv (//)] E. (A) E. (p) 



6><^)[v,£:^](z) dz 
iz -A){z- Id) (liny 



, for A,p € {Au...,An} . (B.43) 



In (IB .431 ) we have used that Ai,. . .,Af^ are all inside of for L large enough. We first expand the N x N 
determinant appearing in the final expression for X]^[v,E'^] into its discreet Fredholm series: 



det 



N 



_ V^'^\Ak,Ai) P^'^\Ak,A() 
Okt + ^^TT-^. — + a- 



1=0 ii, ...,;„ 



det„[v(^)(i,„,i,J + aP(^)(A,„,^,j] 
n\ nLi [^^(^4)] 



(B.44) 



Next, observe that 



{Zin) 



Za - At 



Y\ {4 sin2 [ttv {A J] El {Aj] x det„+i 



a=l 



(Za-Ai^) ^ a 

-iO^^^[y,EZ^]iA„) 1 



(B.45) 



It can be readily seen that for any z belonging to the interior of 

%/^)Cy,y(y)) 



0^'^[v,EZ'](z)= J\ 



2n y - z 



EZ^ (y) with 'ir(^) = '^hk-L U ^T;L U '^^l U % U %d-L , (B.46) 
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and /^^^ is as given by (IB.42I) . Then, using the multilinear structure of a determinant, one gets that 



a 



E-_\y)f%,v{y))det„^i 



-/o(«[v,£:2](^,j 1 

iZa-\)-' 1 



= det„+i 



J In 



n 

EZ\y) f'^Ky, V {y)) Ff det„ 



1 



Za - \ . 

(B. 



This leads to the claim, once upon inserting this representation into the discreet Fredholm series. 
B.5 The Natte series for a Fredholm minor 



In this subsection, we recall the form of the Natte series representation for the Fredholm minor (IB .341 ) involved in 
the representation of form factor sums in generalized free fermionic models. We refer the reader to theorem 2.2 
and proposition 7.2 of reference f6T\ for further details relative to this Natte series expansion. 

Let El = e-'^"('^)-^;('^) be such that 

• u and g are holomorphic in the open neighborhood Us/2 of K.; 

• u has a unique saddle-point Aq on the real axis which is of order 1, ie u"{Ao) < 0; 

• the function y is holomorphic in an open neighborhood Nq c Us/2 of [ ; ^ ]. 

Also, let '4"'^ - ^^'"^ n [z € C : \%z\ < w]. The contours "^^"^ and ^^"'^ have been depicted in Fig.|2l 

For w > \Ao\ + q > and x large enough, the Fredholm minor X,An}[v, E^] defined in (IB. 341 ) admits the below 
Natte series representation 



yv'^(q)+v^{-q) 



2nxu"{AQ) 



+ ■ 



n>l 7C„ £„(^)^^'^' <e-/|*l 



^l+2v(?) 



{linf 



(B.48) 



The +/0"'" regularization of u only matters in the time-like regime (where \Aq\ < q). The functional S, and J?lo 
are given respectively by (14.461) (14.471 ) and (14.481) . The notations and the structure of the sums appearing in the 
second line of (IB.48I ) are exactly as explained in theorem |4~T] 

The Natte series is convergent for x large enough in as much as, for n large enough. 



< C2 



(B.49) 



There c\ and C2 are some ^-independent constants. They only depend on the values taken by u, and g in some 
small neighborhood of the base curve 

•^^^ and by y on a small neighborhood of [ ; ^ ] , whereas 



C3 - - mm 
4 



|l/2, l-2m_ax|5l[v(T^)]|-Y,j where - 2sup | |5l [y(z) - y(T^)]| : k - t^| < 6 , t = ±| . 
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Here e > is sufficiently small but arbitrary otherwise. We stress that, should these norms change, then so would 
change the constants ci,C2 and C3 but the overall structure of the estimates in x would remain. 

The Natte series expansion (IB.48I) has a well defined w +00 limit: all the concerned integrals are convergent 
as the functions //^-j'^ approach zero exponentially fast in respect to any variable that runs to 00 along 
Moreover, this limit does not alter in any way the estimates (IB. 491) ensuring the convergence of the Natte series 
(the constants C1-C3 are w-independent). 

We now list several properties of the functions 

i) {{u (zt)} ; {zt}) [v] is a function of {u (zt)} and {zt). It is also a regular functional of v. 

ii) Hli'f {{u iz,)} ; {z,}) [yv] = O (y") and the O holds in the (L^ n L'°){'^^^^) sense. 

iii) /Z^-^'" can be represented as: 

Hlif {{u (Zt)] ; [zt]) [V] - Hlif ({V (Zt)} ; {u (Zt)} ; {zt}) i>c M (zt))"^" x ]~[ (e~2'-(-) _ 1)' . (b.50) 

with H\l^f is a holomorphic function for |5l (v)| < 1/2 and k [v] {A) is given by (14.471 ). 

iv) One has uff - O (x"") and for « > 2 

b=Q p=Q ,„=h-[^]^ y VXL jy re|±l;0) ^ 

The O (x~°°) appearing above holds in the (L^ nL~)(^^~^) sense. In order to lighten the formula, we have dropped 
the argument-dependent part. However, we do stress that the O {x~°°) as well as [^^1?'^]^ ^ ^ depend on the same 
set of variables as H^^'x- Also, we agree upon t] = 1 for /Iq > ^, = -1 for \Ao\ < q and we made use of the 
shorthand notation 

e+ = e''^"('?)x-2^(«) , e_ = e^^^-^) x^^^-^) and e,, - (1 + »/) e'^"(-^«) . (B.52) 

Finally, the functions [^n;?'^],„y, /,r ^"^^y supported on a small vicinity of the points ±q and Aq. In such a 
case, the contour of integration reduces to an integration for each variable Zt to a small circle dDo-^r around Vr 
(v+ = ±q, vo = Aq). Their dependence on x is as follows. If a variable zt is integrated in a vicinity of v^. The 
function [^l^'^Jmpfor contains a fractional power of x-^'^^^^>^~^'-^'^^\ multiplied by a function of zt which has an 
asymptotic expansion into inverse powers of x. This asymptotic expansion holds on dDo-v^. The coefficients in 
this asymptotic expansions contain poles at zt = Vr- By computing the integrals associated to the terms in this 
asymptotic expansion through the poles at Zt = one obtains that function coefficients associated to x~'' terms 
produce, in fine, a contribution that is a (Inx/x)'^. Finally, the structure of these poles is such that, upon computing 
all the partial derivatives and for any holomorphic function h in the vicinity of the points ±q, Aq, one should make 
the replacement: 

2 eth {zt) ^ Tib Qi (io) - h i-q)) + (m - Tjp) {h (q) - h {-q)) + + +{1+1}) (5,;o/2) h (v^) . (B.53) 

There is one last property which we conjecture to be true for the detailed representation (IB. 511 ) of but 
that has not been proven so far. Namely, 
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Conjecture B.l For a given n the sum in (IB.51I) only contains those combinations of the integers m,p,b and t 
that satisfy to the constraint 

{m-tjp + Sr,if + ^+{m + tj{b- p)- 6r,-i)^ < n . (B.54) 

C Multidimensional Fredholm series for liniA^^+oo p^^^^f (x, 

We begin this appendix by deriving the so-called discreet multidimensional Fredholm series representation for 
P^N-is^-^' ^^^^ prove in theorem IC.ll that this representation has a well defined thermodynamic limit that 

we denote p'^\x, t). This analysis will allow us to provide (proposition IC.ll ) yet another representation for the 
thermodynamic limit p^^{x, t). This alternative representation for p^"^{x, t) is used in subsection 14.61 so as to 
construct the multidimensional Natte series for p^^ix, t). 

Theorem C.l p^!^jj(x, t) admits a well defined thermodynamic limit p|.™?(x, t) that is given by a multidimensional 
Fredholm series 

1 



p^^(,x,t) 2_,c dY''J ilinrfilinrJ ilinf^^ nLAZk-At)iyt-Zk)(yn.i-A,) ^" 



—0 . v^"v . v^^^v . {linf rr,= i 

n n+\ 

n n CVfo - ^ - ic) (/la -Jb- ic) n 

— 1. n -"^i'^" <-«"i) 

n (ja-yb- IC) n {^a - h- IC) 

a,b=l a,b=\ 



Za - h 

det„ \6ke + jVkt det„ [du + yVk[ [yv] ] 



det [/ + jRIln] det^ [/ - Kjln] 



7=0 

(C.l) 



The function f appearing above is supported on the contour ^ = "^^""^ U '^q. The contour is a small loop 
around \ —q',q\ whereas '^q is a small loop around Both "^q and ^q lie below the curve ^^""^ as depicted on 
Fig. |6] All of the aforementioned contours lie inside of the strip Us (12.151 ). The function f is supported on and 
reads 

f(y,y(y))-h^i^>(y)- ^_,J^^^_^ h^^(y) with ^ = ^~)u^. (C.2) 

There 1^ stands for the indicator function of the set A. The function v appearing in the n^^-summand of (IC.ll) 
corresponds to the unique solution of the linear integral equation driven by the resolvent R of the Lieb kernel 
{ie [I - Kiln] [/ + Rlln] = I): 

%, , n n+1 

v{A) + rj ^R{A,p)v{p) = -Z{A)ll + J]0(i,^) - 2<^(i,jJ . (C.3) 

_g £1=1 a=l 

Hence, v depends on the integration variables Aa (with a - I,. . . ,n) and ya (with a = \, . . . ,n + I), ie v{A) 
\'i ! I {Aa}" ; {ja}"''^^ ). We kept this dependence implicit in (IC.ll) so as to shorten the formulae. The entries of the 
finite-size determinants are as defined in (14.81) . They depend on the same set of auxiliary variables as v. Finally, 
we agree upon 

n+l n ?, 

^l({Aa]'l ; [yc^T' I r) = "0 ^y<^^ - Z ^"^-^ + d - r) | «o iA)v{A) dA . (C.4) 

a=l a=l ^„ 



53 



. + i6 



• . 


\ 


R 


-w 




'\ w 








-iS 







Figure 6: The contour consists of the solid line. The contour '^j^-' corresponds to the union of the solid and 
dotted lines. The loop 'r^q is depicted in solid lines whereas the loop is depicted in dotted lines. 



Proof — 

In order to implement the substitution of the operators d^^^ and drjj^, (cf section 1431 ) in the representation (14.411) 
we introduce, exactly as it was done in the proof of theorem 14.11 the functions E- (A) (14.511 ) (whose definition 
involves the functions g - gi^^ + g2,r cf (14.521 )) as well as (14.361 ) and vjr (14.321 ). 

We then consider the discreet Fredholm series representation for ffyX^ [yVi,^?] _^ obtained in proposition 
IB.4I This will allow us to compute the relevant Taylor coefficients {cf subsection |43] equation (14.201 ) and (14.211) ) 
arising in the representation (14.411) for p^'"^„.(x, t). One has that 



fN+l 



a=l 



n=0 



.(-ir 



'lr=o 
d"z 



nl ^ J {liny 



2 f^S^/^^h <«) 



The contours '^'^^\ '^q have been defined in proposition IB. 41 We stress that the summation over n in (IC.5I ) could 
have been stopped atn-m since, prior to taking the y-derivative at y = 0, the n^^ term of the series (IB.41I ) is a 
smooth function of y that behaves as O {y"). We have set 



yn+l - Zj 



,i;rtY<'->(|itl;lwl;l»)l7) 



, (C.6) 



k=i 



the function f ^^{y, v (j)) is given in (IB.42I) and we have set 

n+l n 

- u{iiN+i) + ^ u{Ak)-uQik) 



N 



(C.7) 



k=l 
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Last but not least, 



^^J^|;^^^|i^{ ,,,:':4,,,.,,., ^rbnl(.)r 





W+1 N X 

7=1 y=l ;=1 k«=0 ''l>'=0 



The functional T^^^ is evaluated at the discretization points (c/ definition 14.21 subsection 14.41 and subsection 
14.21 ) tp, p = 1, . . . , 5 for the contour %)ut encircling the compact K2q- %mt has been depicted in the Ihs of Fig. [T] 
The functional F^^^ reads 

N 

T^'^\v}(4D= <f>(M,l^j)-<f>(f^,^j) with /I,- and defined by ^(mj) ^ j/L md ^Aj) ^ j/L . {C.9) 

We do stress that the variables y^, with k = l,...,?i + l, and /ip or Ap with p = I, . . . ,N + I appearing 
in (IC.5I )- (IC.8I ) are all located inside of the compact K2Al, where Al is such that L^(-Ai) = -wl - 1/2. As a 
consequence, the singularities at A = ti^p of the functions ^2,/- (^) (IC.8I ) are always disjoint from the variables y/t,jUp 
or Indeed, t\^p and f2,p with p - 1, . . . , r stand for discretization points of the contour '^outiin appearing in the 
rhs of Fig. [U cf subsection 14.41 These two contours are such that dS^outim^^iAd > uniformly in L. 

According to the prescription that has been adopted in section 14.31 one has to compute the -derivative of 
representation (IC.8I ) prior to implementing the operator substitution. For this, consider any smooth function w (y) 
such that w (y) = O (y") at y = 0. By applying the Faa-di-Bruno formula, we get that 



0=1 ^lr=o \c„\ p=i 



lr=o 



#(i,...,i„ 



n(t) jx^^Kij-^^u- (c-io) 



5^0 

X 



There the / in front of the sums indicates that these are constrained. The first sums runs through all choices of 
N + 2 integers €p>0 such that 

w+i 

= , for p = \, . . . ,n , ij^+i > n and ^ €p = m . (C.ll) 

p=0 

The second sum runs through all the possible choices of sequences of integers kpj with 

p = \,...,N and j = \,...,tp such that j^p,] = ■ (C.12) 

7=1 

Finally, we agree upon 

fp 

\^p\^Yj^P'j aiid have set ^ [/Ip (rp)]|yp=o with defined by ^y^y^{^Ap) ^ p/L . (C.13) 
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In (IC.13I ). we have explicitly insisted on the fact that Ap is a function of the parameter jp. By substituting the 
representation (IC.IOI ) on the level of (IC.5I )- (IC.8I) . one can implement the operator substitution "-^ 5^-^ and 



bj,k ^ dr^jt on the level of (IC.8I ). 

The functional ...,,„ and ^^J^^^ are regular in the sense of definition 14.11 Moreover, as L and hence 2Ai are 

large enough, and defining Up^ is chosen in such a specific wa)ll|that the constant of regularity Cc^ja^ '^^2Al 
satisfies (I4.30I ). one gets that 

{^a\\ ^ n,...J„ [rVs (* I ka))] and {771,,,}'; U {r]2,p}[ ^ ^2A,b^ri* I {lla,p})] (C.14) 

are holomorphic in respect to {^^li £ A^q > Wi.pl'i U|?72,/)li £ ^q"^, where No is an r and s independent neighborhood 
of e C. As the constant of regularity C^ja^ large enough and \y\ can be taken small enough, the size of the 
neighborhood A^o is large enough in order to ensure the convergence of the series of differential operators issuing 
from the exponentials e^'^ and e^2.% once upon the operator substitution is carried out. In virtue of corollary ID. 1[ 
and similarly to the summations (I4.55I )- (I4.56I ). the action of the translation operators can be computed directly 
under the integral sign in (IC.5I ) (the integration contours being Cartesian products of one dimensional compact 
curves) and prior to taking the partial Tp or y-derivatives in (IC.IOI ). There are also the differential operators arising 
from the substitutions Op "-^ d^^^ in (IC.8I ) for those pai^ameters Op that are written down explicitly. The resulting 
(9^^, -derivatives should appear outside of the integrals that are written down in (IC.5I ). However, the integrand of 
these compactly supported integrals is a continuous function of the integration variables that is holomorphic in 
respect to {^pl'j, this uniformly in respect to the integration variables. As a consequence, one can exchange the 
derivation and integration symbols in this case as well. 

Note that the constraints (IC.13I) on the kpj's ensure that in (IC.IOI ) there is at most m-n integers kpj that differ 
from zero. As a consequence, there will be at most m translation operators in respect to the 77, ^ variables to take 
into account once that the operator substitution is made. More precisely, the substitution bj^t d,jjf, shifts the 
parameters Tjj^i^ in rur {A, {rjjk}) (14.321 ) to the below value 

n+l ,11, , N 

\ 1 1 



I ^ I 1^1 1 
T]jk = > > + > 



(C.15) 



p=l p^i r-<p r-JVTi 

where the ulimate sum in (IC.15I ) only involves m terms at most. Under the substitution a,, d^^, the exponentials 



in dC.SI) produce a translation of the function ^ y^ , where 



j=i V j 



a=\ 



After carrying out all these manipulations, we are led to the representation 



in particular it depends on L, cf lemma lAT2l However, \~yP\ ■ L is still very small cf lemma lAT2 
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where -Cw^) is a truncated Lagrange series: 



m] 



I 1 1 



=0 



^ Zj 1 1 1 1 TTI 1 1 IjI7^ 



\p=l 



'lr=o 



0=1 
...,(■„ 



5t 



7;2Ai. 



[ti7,(*|fe))])|,=o ■ (C.18) 



\ra=Ma 



We now take the r ^ +oo hmit of (IC.17I )- (IC.18I ). The very construction of C7,.(/1 | {ijj^k}) along with the choice 



of parameters rji^k given by (IC.15I) associated with the fact that 
large regularity constant, leads to (c/ proof of proposition ID.ll ) 

HI* 



lim [mA = 



is a regular functional with a sufficiently 



(C.19) 



this uniformly in Jq, //q /i^ and Tq belonging to ^2/1l- The function // has been defined in (IA.131I 

This uniform convergence also holds in respect to any finite order partial derivative in these parameters. The 

uniformness of this limit in respect to the integration parameters occurring in (IC.17I) allows one to take it directly 

under the integral sign over a compact domain. 
As a consequence, we get that 



a=\ j=l 



[J'- I 



a=l 



Ta=t^a 



(C.20) 



\ra=o 



To get the rhs of this equality we have, in addition to exchanging the limits and derivatives, applied the Faa-di- 
Bruno formula backwards. The constant of regularity of '^^^2Al being large enough, the action of ^^^2Al ^ 
written in the second line of (IC.20I ) is indeed well defined. 

After collecting the various derivatives into a single one, we arrive to the representation 



(C.21) 



where we have set 



Hi* 



(C.22) 



Since, no confusion is possible on the level of (IC.21I )- (IC.22I ). the y-dependence of the parameters Ap, p - I, . . .,N 
is kept implicit again. We also remind that these are functions oivg. 

The truncated ^-dimensional Lagrange series (IC.21I ) together with its s — > +oo limit has been studied in 
appendix ID. 5. 21 It follows from the latter analysis that the s — > +00 limit is uniform in respect to the parameters 
{{yk}"^^, {Zk}") on which ...,;„ depends. Therefore, this limit can be taken under the integrals signs. Similarly, 
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one can exchange the hmit with the y-derivative symbol. It follows from the results gathered in appendix 
ID.5.2l that 



p=y/''> I |-y=0 



(C.23) 



The answer is expressed with the help of v^^\ the unique solution (for y-small enough) to the non-linear integral 
equation driven by the functional F^^^: 



N 



v(^)U)-O^-l/2)ZU)-0U,^) + J]0U,/i,)-20U,3;,)- Yj '^(^"^") with 



^ = cV)(«/^) 



0=1 



a=\ 



0=1 



(C.24) 



Also, in (IC.23I ). appears the Fredholm determinant of the linear integral operator acting on a small loop around 
[ ; ^ ] whose kernel is given in terms of the functional derivative dT'-^^ [p] Qi) /6p (^). The definition of the 
functional derivative is given in (ID. 111) . 

Lemma IA.4I allows one to reexpress the functional ^^J^^^ appearing (IC.22I ) in the case where the parameters 
Aa and fia are defined exactly as in (IC.24I ) in terms of the unique solution v''^\ 

This leads to the below representation: 



I"). A ,■ V (-1)" a"' V X f"'^ FT I 

') = ft§ ^— f ^ J ^ n 



n=0 ;i, ...,!„ 

i^ell;^]!"* 



{\){zk-\) 



det„ 



1 



Zo - /l/i 



n 



Jn+Y - Zk 



.LrtY<«(|^)f;{A/„)f+';{y„)r'l7) 



n n+ 1 

n n (jfo - - ic) {Ai^ -yb- ic) 

a=l h=l 



det^^ 



(5F(^) [pi 



n+1 



n (3^0 -Jh- ic) n - h, - ic) 

p=yy(L) 0,6=1 a,h=i 



det„ + rHdrv^^^]] det„ [ste + rn-f[ry<^^]]) ({iU';. iyaVr') (C.25) 



|\Cyyt -ZA:)(3'«+l -'ti j 

^ det^+i [H(/^) [^] ] det^v [H('')[fy,,L,]] 

Above we have written down the dependence of both determinants on [AjJ and [ya) as a common argument. 

There is no problem to carry out the analytic continuation in (IC.30I) from e f/gg up to y6 = 0: the potential 
singularities that could appear in the determinants are canceled by the prefactor 0^=1 sin^ |^7ryv^^^(/l4)j. From now 
on, we can thus set /? = 

In order to prove the theorem, it remains to take the thermodynamic limit of (IC.30I ) a.tp = 0. 



L — > +00 behavior of v^^^ 

It was shown in appendix ID. 5. 21 equation (ID.46I) . that v^^^ admits a large L asymptotic expansion v'-^^A) = v(A \ 
{AiJI ; {>'a}"^^ ) + 0(L"^). There the O is holomorphic and uniform in some open neighborhood of the real axis 
and the function v(A) - v{A \ [AiJ" , {jol"^^ ) stands for the unique solution to the hnear integral equation (IC.3I) 
(here we have already set P -0). 

As all of the functions we deal with are smooth functions of v^^\ we are thus able to replace everywhere y^^^ 
by V, up to O (^~') corrections. 
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Building on the large L asymptotics of F^^^ and v^^^ it is shown in subsection ID.5.21 (ID.44I )- (ID.45I) . that 

R 



det<^„ 



= det 



(i+o(z.-)), 



(C.26) 



with a O that has the same uniformness properties as stated before. Above, we did not insist that the Fredholm 
determinant det [/ + yR/ln] conesponds to an action on [ ; ^ ]. 



L ^ +00 limit of 'Z/^^^ 

The thermodynamic limit of K^^\{Aa} ; {//a) ; {ya) I 7) is readily computed by using that 



^ (jUa) - ^yvU-)Ua) ^ ^ -- 



(C.27) 



The remainder O (^~^) is uniform in a e [[ 1 ; A/^ ]] and holomorphic in respect to the variables ya and Za belonging 
to Us/2- By using the Euler-MacLaurin formula, the linear integral equation (IC.3I) satisfied by y and the integral 
representation (12.131) for u one gets that 



(C.28) 



with a O that, again, is uniform and holomorphic in respect to or ya belonging to Us/2- It is also holomorphic 
in > 0. 

By using the densification of the parameters Aa and /i^ on [ ; g ] , it is likewise easy to check that 



detA,+i [E^'> [^] ] det^v [^^''%yi^^]] = det2 [/ - K/2n] • (l + O [l-')) 



(C.29) 



L — > +00 limit of the remaining terms 

It is also readily seen due to the densification of the parameters Aa on [ -<7 ; <? ] that the sums over the discreet 
sets Ai^^ can be replaced by integrals over [-^;^] up to 0(L~') corrections. Finally, it remains to estimate the 
contributions of the functions f'-^\ If one focuses on the contributions of the integrals over ya, a = I, . . . ,n + I 
along the curves ^-f/i;i ^^'^ '^hd;L, then one readily convinces oneself that one deals with the type of integrals 
studied in the proof of proposition IB. II Namely, these are precisely the integrals appearing when deriving the 
estimates for the functional /^^^ given in (IB. 61 ). Clearly, each of these integrals can be estimated successively. 
By repeating word for word the proof given in proposition IB. II one has that each of these integrals produces a 
0(L/wl) = o(l) contribution. Hence, this part of the contour 'i^^^^ does not contribute to the thermodynamic 
limit. 

(m) 



As a consequence one obtains the following representation for pJy.ggC-^, 0- 



(m) 



v^c(-i)" d'" r „ r d"z rd"+v 



det„ 



n=0 
11 n+1 

n n Cvfo - ^ - {Aa -yb- ic) 

a=l b=l 

n (ya-yb- ic) n Ua - Ab - ic) *^=i 

a,b=l a,b=\ 



n {Zk - Ak) iyk - Zk) iyn+i - Ak) 
k=\ 



yn+l - Zb 



f]{4sin2 [7ryv(i^)]) 



det„ [Ske + yVkeiyv]] det„ 



Za - Ab 

Ske + yVkdyy] 



det [/ + yR/2n] det^ [/ - K/2jt] 

x(l+o(l))|y=o ■ (C.30) 
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In order to obtain the representation (IC.ll ) for the thermodynamic hmit it remains to drop the o (1) corrections. 

■ 

We now provide an alternative representation for the thermodynamic limit p^^fix, t). 
Proposition C.l The function Pg"j\x, t) admits the representation 

p^^^{x,t)= Mm hm hm lim:— £2()e-^ X („o yv,, Ucfi t^'] :. (C.31) 

The contour 'lo^^^ = ^^°°''n|z e C : (z)| < wj corresponds to a compact approximation ofrf^^ as depicted 
in Fig.\2\ It is such that Um '^1"'^ = The functional X,aw) has been defined in (IB.34I ). 

w—>+ix c e 6^ 

The functional [cT,.] appearing in (IC.31I ) acfi on f/ie Zoo/3 ^ (^zw) fl^ncf /lai' Z^een defined in lemma lAl] 

The compact approximation "^^"^^ of the contour "^^"^^ appearing in (IC.31I ) is there to ensure the well-definiteness 
of the translation operators. Indeed, in the setting discussed in subsection 14.141 and appendix the translation 
operators are, a priori, only defined for functionals that involve the values of their argument on some compact 
subset of C. As a consequence, a priori, the w +oo limit and r +oo limit do not commute. 

Also, the jS ^ limit and the w +oo limits do not commute. These limit should be understood as follows. 
Given w fixed and large enough, one considers the regular functional '^^2w introduced in lemma IATT] The value 

of w defines an associated /3o £ C and jo > such that '^^2w ^ regular functional for j3 e U/j^ and lyl < yo with a 
regularity constant large enough (in particular satisfying (14.301 )). These % ifio) and yo are such that % (J5q) +oo 
and yo — > when w +oo. 

Proof — 

Let Ez} be as given in (I4.51I )- (I4.52I ). v., as in (14.361 ) and xur (14.321) . In order to implement the operator 
substitution, we first expand the functional X,aw{v, E^] appearing in the rhs of (IC.31I) into a series very similar 
to the one occurring in the proof of proposition IB. 41 The sole exception is that, this time, the sums over Aj^ 's 
are directly replaced by integrals over [-q',(]] of the corresponding variables. Also, the function /*^^^ (resp. its 
associated contour '^'^^^) should be replaced by / (resp. ^^"'^ - ^^^"^ U ^^). At the end of the day, one deals with 
the multi-dimensional Lagrange series below 



,^+c«Zj n\ j (2in)" J (lin)" J (lin)"*^^ ^ ^''^""^ 
The functional T appearing above reads 

"4sin2[;ryv,(/lA:)]" 



^ ^^^^^ ^ n,._,jyy.yr.yn. j-j | y_n^ | . ^^^^^ 
U (zk - W iyk - Zk) 



Za - h 



i.x u\A)yv,(A)dA -ixu(tj) a=l 



n+l 

YY Q-ixu(ya) 

k=y a=l 



And we have set 



[apT"\ d"' \ ngje^ 



X _ /f-' ^TT-e ■ -^[yv.lu-o^^ V^r\\ ■ (C.33) 
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The functional T is evaluated at the discretization points tp, p - 1, . . . , for the contour '^out appearing in the 
Ihs of Fig.d] 

One can implement the operator substitution on the level of (IC.33I ) as it was done in the proof of theorem 
IC.ll The well-foundedness of these manipulations (in particular the justification of the exchange of various limits, 
partial derivatives and integrals over compact contours) is justified along very similar lines. Once upon taking the 
r +00 limit we end-up with the below multidimensional Lagrange series 

=0 

There, we agree upon 



ysU;ka}) = VsU;k«}) + Y ^7:^77— it -^fe'^) + Y^itbUa) - Y(/>(tb,ya) \ ■ (C.35) 



Also, the function cr is to be considered as a functional of 

y-^ - y-^ f ^^dr. 



(C.36) 



The multidimensional Lagrange series (IC.34I ) has been studied in appendix ID. 5. II Its s — > +00 limit is uniform 
in respect to the auxiliary parameters {Za\\ and {ja)"''"^ Hence, just as in the proof of theorem fP. II one is 

allowed to exchange the s — > +00 limit with the integration over the compact contours. One can then apply the 
results of appendix ID.5 . 1 1 leading to 

The function v appearing above is the unique solution to the linear integral equation 

?, , n n+1 

vU) + r I ^/?U,/i)y(//) = 0y3-l/2)ZU) + Y^<l>{A,Xa) - Y^<p{X,ya) . (C.38) 

One can build on this result so as to simplify the obtained expression. The expression for the functional function 
^^^^jTcr[v]] is simplified with the help of lemma IA31 

By using the linear integral equation satisfied by y together with the representation of u in terms of (p and mq 
(12.141 ). we get that the oscillating factor present in T [yv] coincides with the one appearing in theorem fC. It 

n+1 n ?, 

2«(3'«) - Ya""^^^^ - "^^^ - ^ I «'U)vU)d^ = 'W({^)MjjrMr)-2i^Pf . (C.39) 

a=l a=l „ 
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We are thus led to the below representation for the rhs of (IC.31I ) 



yn+l - Zb 



Za - h 



n n CVfo - ^ - ic) (Aa -yh- ic) „ det„ + jVu [yv]] det„ 6u + yVke [yv] 

Y\[4sm'[nyvm] 



. ^ A . 1 , • N 1 1 det [/ + 7/?/2;rl det^ [/ - K/2n] 

a,b=l a,h=\ 

(C.40) 

Tha auxiliary arguments of the entries Vut [v] and Vm [v] are undercurrent by those of v. 

One can carry out the analytic continuation from yS e f/^o up toyS = as the potential singularities of the two 
determinants are canceled by the pre-factors 0^=1 sin^[7r7v(/i<:)]- 

There is no problem to take the w +oo limit of the above integrals. Indeed 'rfj^°°-' is chosen in such a way 
that e"^"^''''\ a = l,...,?i+lis decaying exponentially fast in ya when — > oo along As the rest of the 

integrand is a O (j"), a = l,...,«+lat infinity, the integrals along 'tf^^ are convergent. Once upon taking the 
y6 ^ and the w ^ +oo limits, we recover the representation given in (IC.lb . H 



D Functional Translation operator 

In this appendix, we build a convenient for our purposes representation of a functional translation. Our repre- 
sentation applies to sufficiently regular classes of functionals acting on holomorphic functions. Our construction 
utilizes multidimensional Lagrange series (see eg. yj). 



D.l Lagrange series 



Theorem D.l SB 

Let Do^r = {z e C : \z\ < r). Assume that 

• ^fS^aY]), j = i, - ■ ■ ,s and Jil^aVi) ore holomorphic functions of {^aVi belonging to the Cartesian product 

• there exists a series of radii rj < r such that for = rj, j = I, . . . ,s, one has \(pj ({^a|)| < rj. 
Then, the multidimensional Lagrange series 

1 5"' 



z n 



is convergent and its sum is given by 

f{{Za}) 









dets 


Sjk - T— ¥'y(ka}) 





(D.l) 



Above, {z\,...,Zs) stands for the unique solution to the system zj - <Pj{{Za]) such that \zj\ < r for all j. The 
uniqueness and existence of this solution is part of the conclusion of this theorem. 
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D.2 Some preliminary definitions 

Throughout this appendix, M and K will always stand for two compacts of C such that K c Int (M), M has n 
holes (ie C c M has n bounded connected components) and dM can be reahzed as disjoint union ofn + l smooth 
Jordan curved] y„ : [ ; 1 ] ^ 5M - U;;+}7a ([ ; 1 ]). 
Let /j be a holomorphic function on M and set 

The points ?Q_p correspond to the discretization points for dM associated with the Jordan curves ja, as given in 
definition 14.21 It follows readily that the function A i-^ fs{A \ {ga\) is holomorphic m A e K. Moreover, given any 
holomorphic function v{A,y) e 0" {m x W-^ where Wy is a compact in C'^^' , {y e N, one has that 

fM\{v{ta,p,y)}) r " A^ = v {A, y) + h {A) uniformly in A € ^ and J e W, . (D.3) 

dM 

This convergence holds since (^, A,y) i-> v(^, j) + h (^) / - A) is uniformly continuous on dM X K x Wy. 

We recall that, given a holomorphic function h on M (and hence also on some open neighborhood of M), and 
S a subset of M, we denote \\h\\s = sup^g^ \h {s)\. 

D.3 Pure translations 

We are now in position to establish a representation for translation operators for functional acting on holomorphic 
functions. 

Proposition D.l Let J^[-]{z), z & W, c C^' be a regular functional in respect to the pair {M,K) and let the 
functions fs, v and h as well as the compacts M and K be defined as above. Then, for any [m, ki,..., ki^ € N^^"*"' 



}=i dyj p=\ a=i ' f".p=o )=i dy'^j 



Above, the ■ inside of the argument of indicates the running variable on which the functional !F[-] (y) acts. This 
convergence holds uniformly in (y, j, z) belonging to compact subsets ofDo^yf^ x Int(Wy) X IntCW^), where 

3 ^ Cr nd{dM,K) 

2|M|mxw, + PIIm |5M| + 27rd(5M,/:) ' 

\dM\ stands for the length of dM, d{dM, K)for the distance of K to dM and C<p > is the constant of regularity 
ofT. Finally 2)o,yo = (z e C : |z| < yol- 

Proof — We first consider the case m = and ki = ■ ■ ■ = k( = 0. We assume that s is taken large enough so that 

s n+l 

p=l a=l 



' we remind that satisfies (0) = ja (1) and 7a|[o i [ is injective. 
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Then, for \-y\ < 2yo and kaJ < 2 IIvIImxw,. one has 



(s n+l \ 
sup \ga,p\ + iMm ) X y y VTT — n" - 1^1 ^"P + ii^ii'w 



ta,p+\ -ta,p\ , ,/ I I „,„ \ 2\dM\ 2 

< -Cf- , 



p=i a=i ^" r"- 'a.Pl 



27:d{dM,K) 3 

(D.6) 

Hence, 

(r, ka,p], z)^r [yU* I iz) is holomorphic in (y, z) € Do,2yo x ^o 'Imm ^ ' ^^"^^ 

this for any s large enough. As a consequence, the below multi-dimensional Taylor series is convergent uniformly 
in (y, J, z) € Do,2ro x W'y x and 



n n e'^'"""'^'^- 5- [r/.(* I {Q,,p})] (z)|^^ 
p=i fl=i 

+00 .V n+l (\ (t W"cip cm \ 

Moreover, for any y € Wy an d y e l)o,2ro' one has the bound ||y/v(* I {v(f„,p, j)))||^ + |y| (||v(-, j)||^ + ||/j||^) < 
As a consequence, by (14.141 ) 

\\r[yU* I {vfe,;„j)))](z)-r[yv(*,j) + yM*)](z)||^^^^^,,^^^,^^ 

< roC \\UA I {vfe,;„j))) - v(i,j) -h{A)\\^^^^ ^-^^ , 

due to (ID. 31 ). The norm in the first line is computed in respect to (y, j, z) € i)o,2jo x VKy x W,. The one in the 
second Une in respect to (i, j) € x Wy. We insisted explicitly on the variable-dependence of the functions so as 
to make this fact clear. 

It remains to show that the convergence also holds uniformly on all compacts of "Dq^jq x Int (Wy) x Int (W^) 
when considering partial derivatives in respect to y, ji, . . . ,y[^_ of finite total order. 

One can exchange any such partial derivatives with the Taylor series in (ID. 81 ) in as much as its partial sums 
define a sequence of holomorphic functions that is uniformly convergent on 2)0,270 x x W^. The same arguments 
can be applied to the sequence of holomorphic functions T[yfs{* I |v(fa,p, j)))] (z) . ■ 

Corollary D.l Assume that the conditions and notations of proposition \D.l\ hold. Let ^(^>) = x • • • x '^{^. and 
^(^z) = ^^x - ■ - x^f^be Cartesian products of compact curves in C such that '^(^>) c Int (Wy) and "^"^^^^ c Int (W^). 
Then one has 

+00 s ( ^ fin \ n n cji «+l a,,, 



d'yy d^=z n — T^'^ [rv (*, J) + y/j (*)] (j, z) ■ (D.9) 



f/i/i uniformly in y belonging to compact subsets ofD^^y^. 
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Note that if T depends on a third set of variables belonging to a compact, the results hold as well in respect to 
this third set uniformly on the compact. 

Proof — 

Proposition ID.ll allows one to conclude, in virtue of the uniform convergence of the sequences, that for y 
belonging to compact subsets of Do.yo has the equality 



d'>yld'^zY\—-^r[7v{*,y) + 7h{*)](y,z) . (D.IO) 
J ,_i dy ■ 



The integrals occurring in the first line of (ID. 101 ) are over compact curves and the integrand is smooth in respect 
to the integration variables (y, z) and the auxihary parameters ga,p- As a consequence, the partial fa^p-derivatives 
can be pulled outside of the integration symbols. ■ 



D.4 Weighted translation 

One can generalize the notion of functional translation with the help of multi-dimensional Lagrange series and 
consider more complex objects. For this purpose, we need to introduce some more definitions. Also, from now 
on we only focus on the case of a compact M without holes. 

Let r [•] (ju) be a one parameter family of functionals such that: 

• There exists a constant Cr > such that if v{A,y) is holomorphic in iA,y) e M x Wy, with Wy c and 
ll^ll/fxw,, < Cr then (i, j) i-^ F [v (*, j)] (A) is holomorphic in M x Wy. 

• There exists a contour ^ in Int (K) such that for ||p||^ + ||t||^ < Cp one has 



r\p](ju)-T[T](ju) = J (p-r)(O^|^0u)|v=.d^ + o(||p-T|y . (D.ll) 

5Y [y] (//) j5y (^) will be called the functional derivative of F. This functional derivative is such that, for any 
T holomorphic on M with ||t||^ < Cp, there exists an open neighborhood (^) of the contour appearing 
in (ID. Ill) such that 

(ju, ^ (m) \y=r is holomorphic in (u, ^ € M X ^ (^) . (D. 12) 

• There exists a constant C'^ > such that for ||t||^ + ||y||^ < Cr one has 

\\r[v]m\M<C'rMK and ||F [r] (ju) - F [y] < C;: ||v - t||^ . (D.13) 

The properties of the functional F [•] (A) ensure the solvability of an associated integral equation 
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Lemma D.l Let the compacts M, K and the one parameter family of functional T [•] {A) be as defined above. Let 
h e ^ (M) and r, yo be such that 

I 2\dM\ \ rmin(l,Cr) , / 2\dM\ \ min(l,Cr) 

2yo(r + IWlM) 1 + ^ / r - 2C-yo 1+ < , (D.14) 

\ 2nd{dM,K)l 2(r + C'^) \ 2nd{dM,K)j 2 

Then for \y\ < jq, there exists a unique solution p to the equation p{A) - r[yp{*) + ■yh{*)]{A). This solution is 
holomorphic in (A, y) e M x Do,yo ^f^d such that WpW^ < f- 

Proof — 

Suppose that p and p' are two solutions. Then for lyl < 2yo one has, by construction of yo, that |y| \\p + h\\f^ + 
lyl lip' + h\\f^ < Cr- As a consequence, 

\\P - p'Wm ^ l|r [yip + h)]-T [y ip' + h)]\\^ < 2Cfyo ||p - p% < \\p - p% . (D.15) 

Therefore, p = p' on M, this uniformly in |y| < 2y(). 

In order to prove the existence, one considers the sequence of holomorphic functions on M: ho = h and, for 
n > I hn(A) = h(A) + T [yhn-i {■)]iA). It is readily seen by straightforward induction that, for all « e N and |y| < yo. 

To WKIIk ^ Cr/2 and \\h„+i - /i„||m < \\hn - /j„-iIIm /2 ■ (D.16) 

Hence h„ is a Cauchy sequence in the space of holomorphic functions on Int (M) x Do-yg- It is thus convergent to 
some holomorphic function h on Int(M) x i)o;jo- Since h (A) = h{A) + T[yh] (A), it can be analytically continued 
to a holomorphic function on M x Do;7o- Then, the function p = h - h solves p{A) = T[y{p + h)] (A). It also 
follows that then Upll^v/ < r H 

Proposition D.2 Let f^ be as in (ID.2I ) and assume that the functional F \p] (p.) satisfies to the assumptions given 
above. Let 'F[-](z), with z € VK^ c C^-, be a regular functional in the sense of definition \4.1\ Set 



£r(7,z) = ^lirn^ : Y\e^bM*\^^p^)P')d,, r[yf,(* \ {^p])]{z) 



s—^+oo 

r=l 



: , (D.17) 

f„=0 



where : ■ : indicates that the expression is ordered in such a way that all the partial derivatives appear to the left 
(cf. subsection \4.3i and t^ are the discretization points ofdM. 

Then, there exists yo > such that Xr (t. z) defines a holomorphic function of(y, z) e ^o,yo x Int (W^). 
The convergence of the rhs of (ID. 171 ) to £.r{y,z) is uniform on compact subsets of Do,yo X Int(Wj), and this in 
respect to any partial y or z-derivative of finite order Xr (7» z) is given by 

T [ypl (z) 

-Cr {y, z) = 1^ IFTTTi with p being the unique solution to p{A) = h {A)+Y [yp] {A) . (D. 18) 



det<^ 



/ - tt— — : ('^) 

5v{p) 



\v=yp 



In the denominator appears the Fredholm determinant of the linear integral operator acting on the contour ^ 
with an integral kernel 6T [y] {A) /6v {p)\y^yp. The contour is defined in (ID. 111) . 

Proof — 
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We have, by definition, £r iy, z) - liniv^+oo £s (y, z) with 

(D.19) 



£Ay,z)= 2] flj;^^}- fl r- [r/. (* I M)] a.) • r [r/. (* I M)] U) 

6N 



&=0 



The above series representation for Xi (y, z) corresponds to a particular case of a muhidimensional Lagrange 
series. 

We start by checking the convergence conditions. Let C denote a common constant of regularity for the 
functional T and Y, ie for any v{A,y) e G(^M x IVy) , with Wy c such that llvll^xWy ^ C one has 

!r[v(*,j)](z)€ ^(^3, X W.) and r[v{*,y)]{A) e G [m x W,) . (D.20) 

Then let r > and yo > be as given by (ID. 141 ) but with Cr being replaced with C. Let s be large enough so that 
Za=i k« - ^fl+il ^ 2 \dM\ and |y| < 2yo. It then follows from (ID31 ) that ||y/, (• | < C for l^^l < r. It is also 

easy to see that for l^-^] < r and for any tp e dM, one has |r[y/i.(* | {^p})] {tk)\ < rjl. Therefore, 

• r [y/s.(* I [qp])] {tk), k=l,...,swdT [y/s(* I [qp])] (z) are holomorphic functions of {^,,1 in D^^; 

• for IqI = 3r/4 with A: € U ; ^ 3 one has |r[y/,(* | {^^1)] {tk)\ < rjl < ?>rlA. 

Hence, according to theorem fP.li the multidimensional Lagrange series is convergent and its sum is given by 

£s (y, z) = ^ \ — ^ '-^ ^ (D.21) 



dets 



Sjk-j^T[yM*\{gp})]itj) 



ki,=rp 



where (ti, . . . , t^) is the unique solution to the system tj = r[y/5(* | {Tp})]{tj) with \Tj\ < r for all / 

It is easy to see that, in fact, Xi (y, z) is a uniform limit of holomorphic functions of (y, z) e ^o,2yo ^ ^z- 
Therefore, Xi (y, z) is holomorphic on all compact subsets of Do.iyo ^z- Moreover, there one can permute any 
partial y or z-derivatives with the summations in (ID. 191 ). It is also clear for the previously obtained bounds that, 
-Ci (y, z) is well defined for any s large enough and this independently of the choice of the points t^. 

We now show that its s — > +oo limit exists and then we will compute it. It is readily inferred from the integral 
representation 

Tj= r ^-^ dY (D.22) 

\d=r n fc-r[y/.(*|{^«})]M) 

that Tj = Tj{y), j = I,. . . ,s, solving the system tj = r[y/s(* | {Tp])]{tj), is a holomorphic function of y for 
|y| < yo- Hence, the function ps{A;y) - r[y/v(* | {Tp})] (/I) is holomorphic in {A,y) e M x £)o,yo- Also, by 
construction, 

pAtj;y) = Tj{y) and MMxDo.y, < ■ (^-23) 

Now let p be the unique solution to p(A) = T [yp + yh] (A) with ||p||^ < r, as follows from lemma IdTT] 
Then, keeping the y dependence implicit, we consider 

p{A)-p, (A) = F [y(p + h)] U)-r[y/,(* | {p{tp)})] (i)+r[y/,(* | {p{tp)})] (A) - r[y/,(* | {p,(fp)})] (A) . (D.24) 
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As llpill^ < r, it follows that 

„ „ 2C'ro|5M| 1 

II<A.IIm < cfro ||/,(* I {pitp)}) - M* I {ps{tp)})\\^ < 2 K) "'^ ~ '^■^"^ ^ 2 ~ ^'^"'^ ■ 

Hence, ||p - - iAvIIm - IIP ~ PsWm /2- On the other hand, it follows from (ID.24I ) that 

lip - P,s- - -A-sIIm = ||r [yp + yh]-T [yU* \ {p(tp)]]]\\^ < roCf \\p + h- U* \ {p{tp)})\\^ (D.26) 

Therefore p., converges uniformly to p on M. Hence, in virtue of the regularity of 'F, 

r[jU* I {pAtpM (Z) r [yp] iz) uniformly in (y, z) e ^o.yo x . (D.27) 

s—>+oo 

It remains to compute the limit of the determinant. It follows from the functional derivative property (ID.l II ) 



that 



By expanding the determinant appearing in (ID.21I ) into its discreet Fredholm series we get 
detc 



(D.28) 



v=0 



^jk - /-r[r/v(* I kpmtj)] ^ ^ , = J] ^^ r detp [a,{m„pc)] 



(D.29) 



with 



As[nq,pe) = ^ 



tk+\ - h 



^ 2in{tk-fie) Sviji^) 



r[v + /.(*|{p.M)](o) 



v=0 



d^ 1 5T[v + yp] 



dM 



v=0 



6r [v + yp] 



(Me) 



(D.30) 



v=0 



The above convergence is uniform in (jiq,fi() e x 'rf. Therefore, by elementary estimates, we obtain that 
the determinant of interest does indeed converge to the Fredholm determinant given in (ID.IBI ). this uniformly in 

Irl < ro- _ 

Therefore, we obtain that Xs (y, z) is a sequence of holomorphic functions on Do,yo x Int (W^). that converges 
uniformly. As a consequence, Xr (t> z) is holomorphic on every compact subset of Do,yo ^ I^it (^z) can 
permute any partial-y or z derivative of finite order with the s +00 limit on these compacts. 



D.5 Examples 

We now treat two examples that are of direct interest for the resummation of the form factor series. In the below 
examples, (f> (A, p.) refers to the dressed phase (12.101 ). We remind that it is holomorphic on U^xUs- In the following, 
the compacts K and M <zUs we. such that [ ; ^ ] is contained in their interior. We will also consider functions h 
that are holomorphic on M. 
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D.5.1 r [p] (ju) as a linear functional of p 

Let r [p] (A) = f'^ d^cp (yu, X) p {^) dfi. Then, given a regular functional 'F[-]{z), z e VK^ c C'^s there exists 70 > 
such that for (y, z) e 2)o,yo x 

T\p] iz) 



■Cr (r, z) 



with p (A) 



1 

-J 



dA4>{M,A)p(ju)diu ^h{A) . 



(D.31) 



The limit defining Xr (y> Z) as in (ID.18I) is uniform in respect to such parameters z and I7I < yo- 
Proof — 

In order to apply proposition ID.2I one should check the assumptions on the functional T. It is readily seen 
that, independently of the norm of y and p 

\\r[v]\\M<2q\\dMMxMMK and ||r[v]-r[p]||^<2^||5„^||^^^||v-p||^ and ^ 0") ^ (5^^) (^,^) . 

The validity of the holomorphicity conditions is readily checked by standard derivation under the integral theo- 
rems. One is thus in position to apply proposition ID.2l and the claim follows. ■ 

D.5.2 Non-linear functional F^^^ [p] (ju) 

We now treat the case of the non-hnear functional below 



co) 



^ (oj) - j/L (oj) - j/L) lin 



(D.32) 



There 7 = H 1 ; A^]] \ {/i, . . . , /„), ^ is given by = ^ + F/L and < j/L < D with N/L D. Finally, % is 

a small counterclockwise Jordan curve around [ -<7 ; <7 ] such that Int (A") D 'tfq. Note that /Iq, resp. Aa, appearing 
in (ID. 321) stand for the unique solutions to ^(jia) = a/L, resp. to (Aa) = a/L. 

Proposition D.3 Let J^[-]{z), z € cz C^^ be a regular functional and assume that N, L are large (and such that 
N/L D). Then, there exists yo > such that for (y, z) e 2)o,yo ^ '^^^ ^ large enough 

r[yp(^'](z) 



(5v (0 



(D.33) 



v=yp' 



where p(^) is the unique solution to p^^\A) - h{A) + F^^^[yp'^^-'] (A). This solution is such that 
p^'''>iA)=p{A) + 0{L-'), 



(D.34) 



where p solves the linear integral equation ^/ + ■ p - h and the 0(L ^) is a holomorphic function ofy and A. 
Moreover, this estimate and holds uniformly in |y| < yo and A e Us- Finally, 



6r^ 



v=yT 



1 



linL 4-i ^-yr (0 - j/L 



Above I + R/ln stands for the resolvent of the Lieb kernel acting on [-q',q]- And one has 



det-g; 



1-7 



( 5F<^) [v] 



0") 



det[_,;,] [/ + y/?/27r].|l + o|i)) 



(D.35) 



(D.36) 
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Proof — 

In order to apply proposition ID.2[ we ought to check that F^^^ satisfies to all the necessary conditions. For this 
we observe that 

e ^ inf |^(w) - ^1 > where the inf is taken for oje'iS'q, i€[0;D]. (D.37) 

Then, we choose a constant C > and consider L large enough so that C < eL/2. It then follows that the functions 
(<^) - ill-, for 7 = 1 , . . . , A'^, have no zeroes on 'loq and some immediate neighborhood thereof provided that 
v{z,y) € i?{M X Wy^ with Wy c C'^' and HvH^xw,. < C. It then follows by the derivation under the integral sign 
theorems that F*^^ [v(*, j)] (jj) is holomorphic in (yu, j) e M x Wy. 

In order to establish bounds on F^^^ [v(*, j)] (ji) for v holomorphic and such that ||v||;vfxWj < C, it is convenient 
to represent 

v(a)) 

Ci^) rdt 



+ r ^ ^^'"^ , . (D.38) 



^(oj) -j/L ^yiaj)-j/L L(^(w)-j/L) J L {^,{aj)-j/L) 

As 'i^^ c Int (^), there exists a constant ci > such that for any function y holomorphic on K, one has ||v'||^^ < 
ci ||y||;f. Also, 

inf 1^, (oj) - 5| > 6/2 where the inf is taken for oje'rfg, s e [ ; D ] and \t\ < eL/2 . (D.39) 
Hence, for any v € ff(^Mx W^) such that ||v||/s:xw_,, < C < eL/2 

MkxW, ' 



ell? 



< II0IImxm^;[^ci|i + ^[||^||j, + 6/2]| IMI^xw,, . (D.40) 
This provides an estimate for the constant C^^^, entering in the bounds for ||F^^^ [^^IImxiv,' ^^^'^ 



-(i)r 

VJ V^^ ~ ^ L' J VA'/ — 

Where, 

p(w) 



this uniformly in /i e M and L large enough. Therefore, 



(^,(c^)-j/L)3 

t(w) 



^ -0(|b-T||y , (D.42) 



It follows that there exists a sufficiently small open neighborhood {^^^ of such that the functional derivative 
is holomorphic in (//, ^) e M x "V (^q)- Moreover, we get that there exists an L-independent constant C2 such that 
||F(^)M-F(^)[r]||^<C2||p-r||j, 
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We are now in position to apply proposition ID.2I It follows that Xp^) can be expressed in terms of the unique 
solution p*^^ to p^^^Qj.) = h{ii) + Y^^^yp^^^] (ji) with ||p''^''||^ < r uniformly in lyl < jq. 
This means that, 



5v{0 



y(d^<p){H,0 j uniformly in (/i,Oe^,'- (D.44) 



v=p'-' 

-q 



In this limit, the contour integral in the Fredholm determinant can be computed and since the Fredholm deter- 
minant of a trace class operator is continuous in respect to the trace class norm (which is bounded by the sup norm 
in the case of integral operators acting on compact contours) 

det^, [/ - 5r<^)[ry] ill) IdviO |,=pa)] ^-^^ det^.^.^j [/ - y^A<p{^i,A)\ = det[_,^.^] [/ + yR/2n] . (D.45) 

Where R is the resolvent of the Lieb kernel. 

We now characterize the leading behavior of the solution p^^'' when N,L ^ +00. By repeating the type of 
manipulations canied our previously, and using that p'-^^ is bounded on K uniformly in L, we get that the non- 
linear integral equation for p^^^ takes the form 



V) — + o - 



There the O is uniform in p e Us- The Riemann sum can be estimated by using the Euler-McLaurin formula and 
the uniform boundedness of p*^^ on K. After carrying out the resulting contour integral over 'loq we obtain 

1 

p(^\/i) = h{p.) - y ^ ^Ripi, s)p^'^\s) + O . (D.46) 



The O appearing in (ID.46I ) is holomorphic in ju e Ug. Indeed, p^^^ just as all the other terms in (ID.46I ) are 
holomorphic on Us- This proves that p*^^ admits an asymptotic expansion in L such that p^^\cjj) = p (a>) -1- O 
where p is the solution to the integral equation (/ -1- yR/2n) ■ p = h . As p'-^^ and p are both holomorphic in Us, so 
is the remainder. Moreover, one can convince oneself that this O is uniform in ju € Us- Therefore, the regularity 
of the functional T [•] (z) (I4l4l) leads to Tiyp^^''] (z) - T [yp] (z) + o(l-^)- m 
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